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Abstract 

We classify all spherically symmetric perfect fluid solutions of Einstein's 
equations with equation of state p = a^i which are self-similar in the sense 
that all dimensionless variables depend only upon z = r/t. For a given 
value of a, such solutions are described by two parameters and they can 
be classified in terms of their behaviour at large and small distances from 
the origin; this usually corresponds to large and small values of z but (due 
to a coordinate anomaly) it may also correspond to finite z. We base our 
analysis on the demonstration (given elsewhere) that all similarity solu- 
tions must be asymptotic to solutions which depend on either powers of z 
or powers of Inz. We show that there are only three similarity solutions 
which have an exact power-law dependence on z: the flat Friedmann solu- 
tion, a static solution and a Kantowski-Sachs solution (although the latter 
is probably only physical for a < 0). For a > 1/5, there are also two fam- 
ilies of solutions which are asymptotically (but not exactly) Minkowski: 
the first is asymptotically Minkowski as z — > oo and is described by one 
parameter; the second is asymptotically Minkowski at a finite value of z 
and is described by two parameters. 

As a simple application of this classification, we then present a complete 
analysis of the dust (a = 0) solutions, since these can be written down 



explicitly and elucidate the link between the z > and z < solutions. 
As in the a/fl case, the most general dust solution is described by two 
parameters. The first one (E) corresponds to the asymptotic energy at large 
\z\, while the second one (D) specifies the value of z at the singularity which 
characterizes such models. The 1-parameter family with z > and D = 
are inhomogeneous cosmological models which expand from a Big Bang 
singularity at t = and are asymptotically Friedmann at large z; models 
with E > are underdense and expand faster than Friedmann, while those 
with E < recollapse to black holes and contain another singularity. The 
D = solutions with z < are just the time reverse of the z > ones. The 
2-parameter solutions with D > again represent inhomogeneous models 
but they necessarily involve both z < and z > regimes, the Big Bang 
singularity is at \z\ = 1/D and (while there is no exact static solution 
in the dust case) they are asymptotically "quasi-static" at large \z\. The 
solutions with E > expand or contract monotonically, whereas the ones 
with E < recollapse to or expand from a second singularity. Depending 
on the values of E and D, this may be either a black hole singularity or 
a naked singularity. The D < models either collapse to a shell-crossing 
singularity and become unphysical or expand from such a state. 

We then discuss solutions with pressure. These share many of the 
characteristics of the dust solutions but they also exhibit new features. At 
large distances from the origin, we show that there is a 1-parameter fam- 
ily of solutions which are asymptotically Friedmann, a 1-parameter family 
of solutions which are asymptotically Kantowski-Sachs and a 2-parameter 
family which are asymptotically quasi-static. All these solutions can be de- 
scribed by parameters E and D, analogous to those which arose in the dust 
case. The asymptotically Minkowski solutions, which arise for a > 1/5, are 
discussed in detail elsewhere. The possible behaviours at small distances 
from the origin depend upon whether or not the solutions pass through a 
sonic point. If the solutions remain supersonic everywhere, the origin cor- 
responds to either a black hole singularity or a naked singularity at finite 
z, as in the dust case. However, if the solutions pass into the subsonic re- 
gion, they reach z = and their form is restricted by the requirement that 
they be regular at the sonic point. There is again a 1-parameter family of 
asymptotic Friedmann solutions: this includes a continuum of underdense 
solutions and discrete bands of overdense ones; the latter are all nearly 
static close to the sonic point and exhibit oscillations. There is also a 1- 
parameter family of asymptotically Kantowski-Sachs solutions. However, 
there are no asymptotically static solutions besides the exact static solution 
itself. 
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1 Introduction 



Self-similar models have proved very useful in General Relativity because the 
similarity assumption reduces the complexity of the partial differential equations. 
Even greater simplification is achieved if one has spherical symmetry (Cahill & 
Taub 1971) since the governing equations then reduce to comparatively simple 
ordinary differential equations. In this case, the solutions can be put into a form 
in which every dimensionless variable is a function of some dimensionless combi- 
nation of the cosmic time coordinate t and the comoving radial coordinate r. In 
the simplest situation, a similarity solution is invariant under the transformation 
r — > at,t — > at for any constant a. Geometrically this corresponds to the existence 
of a homothetic Killing vector and is sometimes termed similarity of the "first" 
kind. We confine attention to such solutions in this paper. We shall also focus 
on the case in which the source of the gravitational field is a perfect fluid with 
an equation of state of the form p = a/j,. Indeed, Cahill & Taub (1971) showed 
that this is the only barotropic equation of state compatible with the similarity 
assumption. We will assume \a\ < 1, as required by causality, and usually require 
a to be positive. Note that "geometric" self-similarity (a property of the metric) 
and "physical" self-similarity (a property of the fluid) coincide for a perfect fluid 
but this need not be the case in general (Coley & Tupper 1989). 

What makes such solutions of more than mathematical interest is the fact 
that they are often relevant to the real world. For example, an explosion in a 
homogeneous background produces fluctuations which may be very complicated 
initially but which tend to be described more and more closely by a spherically 
symmetric similarity solution as time evolves (Sedov 1967). This applies even if 
the explosion occurs in an expanding cosmological background (Schwartz et al. 
1975, Ikeuchi et al. 1983). The evolution of cosmic voids may also be described 
by a similarity solution at late times (Bertschinger 1985). The same idea applies 
in a wide range of contexts in fluid dynamics since, in this case, self-similar 
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asymptotics can be obtained from dimensional considerations. Indeed it was in 
this context that the concept of self-similarity of the first kind was first introduced 
(Barenblatt & Zeldovich 1972). Recently it has become clear that spherically 
symmetric self-similar solutions also play a crucial role in the context of the 
"critical" phenomena discovered in gravitational collapse calculations (Choptuik 
1993, Evans and Coleman 1994, Gundlach 1995, Koike et al. 1995, Maison 1996). 

In the cosmological context these considerations led Carr (1995) to propose 
the "similarity hypothesis". This says that - under certain circumstances (eg. 
non-zero pressure and high non-linearity) - cosmological solutions may naturally 
evolve to a self-similar form even if they start out more complicated. There 
is evidence for this in both the spherically symmetric context (Carr & Coley 
1998a) and the spatially homogenous context (Wainwright & Ellis 1997). The 
possibility that self-similar models may be singled out in this way from more 
general spherically symmetric solutions means that it is essential to understand 
the full family of such solutions. Unless the pressure is zero, similarity solutions 
generally have a shock or sonic point and the nature of the solution at this point 
plays a crucial role. The original study of Cahill and Taub (1971) focussed on 
solutions with shocks, whereas subsequent authors (eg. Bogoyavlenski 1977) have 
focussed on solutions with sound-wav! es. In this paper we will only consider the 
latter and we will focus exclusively on solutions which are "regular" at the sonic 
point in the sense that they have a finite pressure gradient and can be continued 
beyond there. Even some of these solutions will turn out to be unphysical. 

Due to the existence of several preferred geometric structures in self-similar 
spherically symmetric models, a number of natural approaches (i.e. coordinate 
systems) may be used in studying them (Bogoyavlensky 1985). The three most 
common ones are the "comoving" , "homothetic" and "Schwarzschild" approaches. 
In the comoving approach, pioneered by Cahill & Taub (1973) and employed by 
Carr and Henriksen and coworkers, the coordinates are adapted to the fluid 4- 
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velocity vector. This probably affords the best physical insights and is the most 
convenient one with which to study the solutions explicitly. In the homothetic 
approach, used by Bogoyavlensky and coworkers, and adopted more recently by 
Brady (1994) and Goliath et al. (1998a, 1998b), the coordinates are adapted to 
the homothetic vector. In this case, the governing equations reduce to those of an 
autonomous system and so dynamical systems theory can be exploited to study 
the equations mathematically. The "Schwarzschild" approach, adopted by Ori & 
Piran (1990) and Maison (1996), is useful if one wishes to match a self-similar 
interior region to a non-self-similar asymptotically flat exterior region. This is 
because one can analyse null geodesies most simply in these coordinates, enabling 
the causal structure of spacetime to be studied. The relationship between these 
different approaches is discussed in more detail in Appendix A. All of them are 
complementary and which is most suitable depends on what type of problem one 
is studying. In this paper it is most convenient to use the comoving approach. 

The intention of this paper is to provide a complete classification of perfect 
fluid spherically symmetric self-similar solutions. We achieve this by analysing 
all possible behaviours at large and small distances from the origin. In the sim- 
plest situation this just corresponds to large and small values of the similarity 
variable z but the analysis is complicated by the fact that (due to a coordinate 
anomaly) a finite value of z may sometimes correspond to zero or infinite distance 
from the origin. For this reason some of the similarity solutions we discuss were 
missed in previous treatments (and indeed in an earlier version of this paper). A 
more rigorous demonstration that our classification is complete is given elsewhere 
(Carr & Coley 1998b) and consists of two parts: (1) an analysis of all solutions 
whose asymptotic behaviour is associated with large or small values of \z\ and a 
demonstration that these always have a power-law dependence on z; (2) an anal- 
ysis of solutions whose asymptotic behaviour is associated with a finite value of 
z and a demonstration that these have a power-law dependence on Inz. We will 
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use this "power-law" property as the starting point of the present analysis. This 
simplifies the discussion considerably and allows us to focus on the nature and 
physical significance of the solutions. Some of them have been found before [see 
Coley (1997), Carr (1998) and Carr & Coley (1998a) for recent reviews], so this 
will require some overview of previous work. However, this is the first time they 
have all been brought together, with the connection between them being made 
explicit. This work complements the analysis of Goliath et al. (1998a, 1998b), 
which also delineates the different types of solutions but without making their 
physical significance clear. The precise relationship between our two approaches 
is discussed in a separate paper (Carr et al. 1998). 

We will show that perfect fluid self-similar spherically symmetric solutions 
have four possible behaviours at large distances from the origin. They are either 
asymptotically Friedmann, asymptotically "quasi-static" , asymptotically Kantowski- 
Sachs or asymptotically Minkowski, with the last family being subdivided into 
two (one of which is associated with a finite value of z). The possible behaviours 
at small distances depend upon whether or not the solutions pass through a sonic 
point and reach z = 0. If the solutions remain supersonic everywhere, the ori- 
gin is at finite z and corresponds to either a black hole singularity or a naked 
singularity; in either case, the small-scale behaviour is uniquely determined by 
the large-scale behaviour. If the solutions pass through a sonic point, they may 
be discontinuous there and the situation is more complicated. However, in this 
paper we confine attention to solutions which are regular at the sonic point and 
physically realistic throughout the subsonic regime; this excludes solutions which 
have shocks or enter a negative mass regime. All such solutions reach z = and 
are either asymptotically Friedmann, exactly static or asymptotically Kantowski- 
Sachs for small z. Which of these situations applies depends on the value of z at 
the sonic point and this itself is determined by the large-scale behaviour. In all 
cases, one can therefore classify solutions by their form at large distances from 
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the origin. By way of introduction, we now briefly describe these forms. 

The first class of solutions is a 1-parameter family asymptotic to the flat 
Friedmann solution at large values of \z\. Attention originally focussed on models 
containing black holes because there was interest in whether black holes could 
grow at the same rate as the particle horizon. Carr & Hawking (1974) showed 
that such solutions exist for radiation (a = 1/3) and dust (a = 0) but only if 
the universe is asymptotically rather than exactly Friedmann (i.e. there is no 
solution in which a black hole interior is attached to an exact Friedmann exterior 
via a sound-wave) and this has the important implication that black holes formed 
through purely local processes cannot grow as fast as the Universe. [In fact, all 
subsonic solutions which can be attached to an exact Friedmann model via a 
sound-wave are unphysical in the sense that, as one goes inward from the sonic 
point, they either enter a negative mass regime or reach another sonic point at 
which the pressure gradient diverges (Bicknell & Henriksen 1978a).] Carr (1976) 
and Bicknell & Henriksen 1978a) then extended this result to a general < a < 1 
fluid, while Lin et al. (1976) and Bicknell & Henriksen (1978b) considered the 
case of a stiff fluid (a = 1). It seems likely that all black hole solutions are 
supersonic everywhere, though this has still not been rigorously proved. 

There are also perturbations of flat Friedmann models which represent density 
fluctuations growing at the same rate as the Universe (Carr & Yahil 1990). These 
solutions are asymptotically Friedmann at both large and small values of \z\ and 
regular at the sonic point. Such transonic solutions can be either underdense or 
overdense relative to the exact Friedmann model. The underdense solutions may 
be relevant to the existence of large-scale cosmic voids (Carr & Whinnett 1998), 
while the overdense ones may be relevant to the formation of highly overdense 
blobs at a cosmic phase transition. While there is a continuum of underdense 
solutions which are regular at the sonic point, regular overdense solutions only 
occur in successive and very narrow bands. These solutions have the characteristic 
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that they are all approximately static near the sonic point, although they depart 
from the static solution as they approach the origin. 

The second class of models is associated with the Kantowski-Sachs solution. 
This is a type of homogeneous model first studied by Kantowski and Sachs (1966) 
for the a = case and by Collins (1977) for arbitrary a. For each a there is a 
unique self-similar Kantowski-Sachs solution and there also exists a 1-parameter 
family of solutions asymptotic to this at both large and small values of \z\ (Carr 
& Koutras 1992). Solutions with — 1/3 < a < 1 are probably unphysical because 
they are tachyonic and the mass is negative. Solutions with —1 < a < —1/3 
may be more physical since they avoid these features. Although such equations 
of state violate the strong energy condition, they could could well arise in the 
early Universe due to inflation or particle production effects. Such models may 
be related to the growth of p > bubbles formed at a phase transition in a p < 
cosmological background (Wesson 1986, 1989, Ponce de Leon 1988, 1990). Note 
that this is the only context in which we will consider negative values of a. 

The third class of models are asympototic to either a self-similar static model 
or what we term "quasi-static" models at large values of \z\. There is just one ex- 
actly static self-similar solution for each (positive) value of a (Misner & Zapolsky 
1964) and there is a 1-parameter family of solutions asymptotic to this. However, 
we will show that there is a 2-parameter family of solutions which are asymp- 
totically "quasi-static" in the sense that they have an isothermal density profile 
at large values of \z\. Since the asymptotically Friedmann and asymptotically 
Kantowski-Sachs solutions are described by only one parameter, such solutions 
play a crucial role in understanding the full family of similarity solutions. Some of 
these solutions remain supersonic everywhere and reach a singularity at finite z. 
However, unlike the Friedmann case, the singularity may be naked. The solutions 
which reach a sonic point may be attached either to the exactly static solution 
or the asymptotically Friedmann solutions. Some asymptotically quasi-static so- 
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lutions have been studied before (Ori & Piran 1990, Foglizzo & Henriksen 1993). 
In particular, they are known to be associated with the formation of naked sin- 
gularities in spherically symmetric collapse (Henriksen & Patel 1991, Lake 1992, 
Joshi & Dwivedi 1993). However, the precise relationship of these solutions to 
the more general quasi-static family has not been discussed before. 

The fourth class of solutions, which only exist for a > 1/5, are asymptotically 
Minkowski and have not been previously analysed at all. They were originally 
found numerically by Goliath et al. (1998b) and this led us to "predict" them 
analytically. There are actually two such families and they are described in more 
detail elsewhere (Carr et al. 1998). Members of the first family are described by 
one parameter and are asymptotically Minkowski as \z\ — > oo; members of the 
second family are described by two parameters and are asymptotically Minkowski 
as z tends to some finite value (though this corresponds to an infinite physical 
distance). The existence of these solutions has been obscured until now because 
one needs to perform a complicated coordinate transformation in order to put the 
metric into an explicitly Minkowski form. As in the asymptotically Friedmann 
and asymptotically quasi-static cases, these solutions may either be supersonic 
everywhere (in which case they contain a black hole or naked singularity) or they 
may be attached to z = via a sonic point (in which case they are asymptotically 
Friedmann at small z). The latter solutions are particularly interesting because 
they are both associated with the occurence of critical phenomena, as discussed 
in more detail by Carr et al. (1998) and Carr, Henriksen & Levy (1998). 

The plan of this paper is as follows. In Section 2 we will introduce the relevant 
equations and discuss the crucial role of the sonic point. In Section 3 we will 
analyse the possible behaviours at large and small distances from the origin, 
emphasizing the key role played by the power-law and log-power-law solutions. 
In Section 4 we will analyse the dust (p = 0) solutions since the complete family 
of such solutions can be derived analytically. Despite the simplifications entailed 
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in dropping pressure, we will find that many of their features carry over to the 
a 7^ case. In particular, the dust solutions illuminate the connection between 
models with positive and negative z. In Section 5 we will consider the solutions 
with pressure. We will show that some of the features of these solutions can 
be understood by combining the insights gained from the dust solutions in the 
supersonic regime with those gained from the asymptotically Friedmann solutions 
in the subsonic regime. However, other kinds of solution have no analogue in the 
dust case. We make some final remarks in Section 6. 



2 Spherically Symmetric Similarity Solutions 



In the spherically symmetric situation one can introduce a time coordinate t 
such that surfaces of constant t are orthogonal to fluid flow lines and comoving 
coordinates (r, 9, <fi) which are constant along each flow line. The metric can then 
be written in the form 



where z/, A and R are functions of r and t. For a perfect fluid the Einstein 
equations are 



where /j,(r, t) is the energy density, p(r, t) the pressure, £/ M = (e~ u , 0,0,0) is the 
comoving fluid 4- velocity, and we choose units in which c — G — 1. The equations 
have a first integral 



and this can be interpreted as the mass within comoving radius r at time t: 



ds 2 = e 2v dt 2 - e 2X dr 2 - R 2 dtf, dtt 2 = d9 2 + sin 2 9 d<\> 2 



(2.1) 



(2.2) 




(2.3) 




(2.4) 
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Unless p — 0, this quantity decreases with increasing t because of the work done 
by the pressure. One can also express it as 

,dR 



it )=47rj o pR 2 —dt' (2.5) 



and this is the more appropriate expression when there is no spatial origin (as in 
the Kantowski-Sachs solution). Eqn (2.3) can be written as an equation for the 
energy per unit mass of the shell with comoving coordinate r: 

*-§0--i>-^-=. —-(f). -e-(f). (,a) 

This can be interpreted as the sum of the kinetic and potential energies per unit 
mass. Only in the p = case are E and T conserved along fluid flow lines. 

By a spherically symmetric similarity solution we shall mean one in which the 
spacetime admits a homothetic Killing vector £ that satisfies 

+ Zuw = 2 9lu ,. (2.7) 

This means that the solution is unchanged by a transformation of the form t — > at, 
r — > ar for any constant a. Solutions of this sort were first investigated by Cahill 
& Taub (1971), who showed that by a suitable coordinate transformation they 
can be put into a form in which all dimensionless quantities such as u, A, E and 

S=?L, M = ^ P = pR\ W = fiR 2 (2.8) 

are functions only of the dimensionless variable z = r/t. Then we have 

- = ---, f?- = £^-, (2.9) 
dt r dz' dr r dz ) 

so the field equations reduce to a set of ordinary differential equations in z. 
Another important quantity is the function 

V(z) = e x ~ u z, (2.10) 
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which represents the velocity of the surfaces of constant z relative to the fluid. 
These surfaces have the equation r = zt and therefore represent a family of 
spheres moving through the fluid. The spheres contract relative to the fluid for 
z < and expand for z > 0. This is to be distinguished from the velocity of the 
spheres of constant R relative to the fluid: 

This is positive if the fluid is collapsing and negative if it is expanding. Special 
significance is attached to values of z for which |V| = 1, Vr — and \Vr\ = 1. 
The first corresponds to a Cauchy horizon (either a black hole event horizon or a 
cosmo logical particle horizon), the second to a stagnation point, and the third to 
a black hole or cosmological apparent horizon. We show shortly that the existence 
of an apparent horizon is also equivalent to the condition M — 1/2. 

The only barotropic equation of state compatible with the similarity ansatz 
is one of the form p = a/i (—1 < a < 1). It is convenient to introduce a 
dimensionless function x(z) defined by 

x(z) = ^fir 2 )- a l {l+a \ (2.12) 

[Note that the factor of An is omitted in the definition of x given by Carr & 
Yahil (1990) but it is required for consistency with eqn (2.3).] The conservation 
equations T^" = can then be integrated to give 

e" = f3xz 2a/{1+a) (2.13) 

A _ „,-l/a c2 



e = 7X 



1/a S 2 (2.14) 



where (5 and 7 are integration constants. The remaining field equations reduce 
to a set of ordinary differential equations in x and S: 

^ +s+ (r!^-^)i s+ < 1+a ^ - < 2 ' 15 ' 
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f2af\ s , + 2S {2 _ 2a)/az{2 _ 2a)/{1+a) _ ^ 4 x W _ \ (1 + a)a . ( i_ a)/aj 
V 1 I y o X \ ^ / 

(2.16) 



1 + + 



(2.17) 



M 

2 



± + ^- 2 x- 2 z 2 ^^S 2 - l -^x-^S & (l + ^j , (2.18) 
where the velocity function is given by 

V = (/3 7 )- 1 x (1 - a)/Q ^ 2 ^ {1 - Q)/(1+Q) (2.19) 

and an overdot denotes zd/dz. The other velocity function is 



while the energy is 



(2.21) 



2 

Eqn (2.18) can then be written in the form 

M = l+( E + \)( V R- 1 )> (2-22) 

so the condition M — 1/2 implies that |Vr| = 1 (corresponding to an appar- 
ent horizon). Note that M ^ 1/2 in the special case E — —1/2 because this 
corresponds to the Kantowski-Sachs solution and, for this, Vr diverges. 

As discussed by Carr & Yahil (1990), we can best envisage how these equations 
generate solutions by working in the 3-dimensional (x, S, S) space. At any point 
in this space, for a fixed value of a, eqns (2.17) and (2.18) give the value of 
z; eqn (2.16) then gives the value of x unless \V\ = -fa and eqn (2.15) gives 
the value of S. Thus the equations generate a vector field (x, S, S) and this 
specifies an integral curve at each point of the 3-dimensional space. Each curve is 
parametrized by z and represents one particular similarity solution. This shows 
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that, for a given equation of state parameter a, there is a 2-parameter family of 
spherically symmetric similarity solutions. 

In (x, S, S) space the sonic condition V = \/a specifies a 2-dimensional surface 
because eqns (2.17) to (2.19) allow one to express S in terms of x and S. The 
same surface corresponds to the condition V = —\/a. Where a curve intersects 
this surface, eqn (2.16) does not uniquely determine x, so there can be a number 
of different solutions passing through the same point. However, integral curves 
intersect | V | = \fa in a physically reasonable manner only if 

f^f.) S 4 + A| x (2-2aya z (2-2a)/(l + a) = (1 + a)a .(l-«)/« (333) 

since otherwise the value of x and hence the pressure, density and velocity gradi- 
ent diverge there. Since eqn (2.23) corresponds to another 2-dimensional surface 
in (x, S, S) space, this will intersect the surface \V\ = y/a on a line Q. Only inte- 
gral curves which hit the sonic surface on this line are "regular" in the sense that 
they can be extended beyond there. (All other solutions would have to contain 
shock- waves.) From each point on this line there will be regular integral curves 
with decreasing and increasing z. One can join any member of the first kind to 
any member of the second kind to obtain a complete similarity solution. 

Physically reasonable solutions cannot have an arbitrary value of x at |V| = 
yfa. If we require 'x to be finite there, then the equations permit just two values 
of x at each point of the line Q and there will then be two corresponding values 
of V. If the values of x are complex, corresponding to a focal point, then the 
solution will be unphysical. If they are real, at least one of the values of V must 
be positive. If both values of V are positive, corresponding to a nodal point, the 
smaller one is associated with a 1-parameter family of solutions, while the larger 
one is associated with an isolated solution. If one of the values of V is negative, 
corresponding to a saddle point, both values are associated with isolated solutions. 
This behaviour has been analysed in detail by Bogoyavlenski (1977), Bicknell & 
Henriksen (1978), Carr & Yahil (1990) and Ori & Piran (1990). 
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One can show that there is a 1-parameter family of regular solutions (i.e. a 
node) only on a restricted part of the line Q and, in the V(z) diagram, this 
corresponds to two ranges of values for z. One range {z\ < z < z 2 ) lies to the 
left of the Friedmann sonic point zp and includes the static sonic point zs- The 
other goes from some value z 3 to infinity and includes zf- We will argue later 
that any family of regular solutions which is described by just one parameter 
asymptotically probably has to hit the sonic line in these ranges. One has a 
saddle point for z < z 1 and a focal point for z 2 < z < z 3 . These features are 
indicated in Figure (la). The values of zi, z 2 and z% can be given in terms of a 
but the expressions are very complicated, so we do not give them explicitly. The 
ranges for a — 1/3 are indicated in Figure (lb); in this case, z% = zf and z 2 = zs- 
On each side of a nodal sonic point, x may have either of its two possible values. 
If one chooses different values for x, there will be a discontinuity in the pressure 
gradient. If one chooses the same value, there may still be a discontinuity in the 
second derivative of x. Only the isolated solution and a single member of the 
1-parameter family of solutions are "analytic" or at least C°° in the sense that 
derivatives of every order are continuous. This contrasts with the case of a shock 
where x is itself discontinuous. 

3 Asymptotic Behaviour of Similarity Solutions 

The key step in providing a complete classification of spherically symmetric per- 
fect fluid similarity solutions is an analysis of their possible asymptotic behaviours 
and we now present this. For simplicity we will assume z > throughout this 
section but the analysis can be trivially extended to the z < case. We will also 
assume a > except in the Kantowski-Sachs case. The full technicalities of the 
asymptotic analysis are presented elsewhere (Carr & Coley 1998b). For present 
purposes it suffices to note that all similarity solutions depend on powers of z at 
large and small values of \z\ or on powers of Inz at finite z. The last possibility 
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arises because, due to a coordinate anomaly, a finite value of z may sometimes 
correspond to zero or infinite physical distance. In this section we will identify 
these asymptotic states explicitly. We will show that there are three exact power- 
law solutions: the flat Friedmann solution, a Kantowski-Sachs solution and a 
static solution. We will also show that, for a > 1/5, there are asymptotically 
(but not exact) Minkowski solutions which asymptote either to infinite z or finite 
z. Finally there are solutions whose origin corresponds to a singularity at finite 
z. The validity of these results is confirmed by the dynamical systems analyses 
of Bogoyavlensky (1985) and Goliath et al. (1998a, 1998b). In particular, the 
existence of the monotone and Dulac functions found in these analyses forbids the 
existence of periodic orbits and limit cycles and thereby excludes other possible 
asymptotic behaviours. 

3.1 Power-Law Similarity Solutions 

In order to find the asymptotically "power-law" solutions explicitly, we look for 
solutions to the field equations of the form 

x = x Q z a , S = S Q z b (3.1) 

where x Q , S Q , a and b are constants. Note that S/S = b and x/x = a. Eqn (2.15) 
is satisfied if 

_ ba[3(b + 1) + a(3b + 1)] 

(l + a)[l + (l + a)6] • 1 ] 

The factor [1 + (l + a)b] cannot be zero since this would be inconsistent with eqns 

(2.17) and (2.18). Eqn (2.16) can then be written in the form 

Az p + Bz q + C = (3.3) 

where 

A - 6[(a-l) + (l + a)(2a-l)6] 2 (i- a )/a o- 4 R _ n , „ w (i-«)/<» q -4 
A= /F(l + a)[l + (l + a)6] X ° 6 ° ' ^ = "(1+^0 \, 

= a 7 2 (fe + l)[2 + 36(l + a)] 

" (l + a)[l + (l + a)6] 1 ' ' 
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and the exponents are 

„ s2o (I^)-« + 2 (I^), , S0 (I^)-4* ,3,) 

Since B cannot be zero, there are then three ways in which eqn (3.3) can be 
satisfied to leading order as z — > or z — > oo and we discuss these in turn. 



• p = q, A + B = 0. In this case, the condition p = q implies 

2a 



l + a 



and the condition A + B = implies 



x 



Eqn (3.2) then requires 



(l-a)/a = io2 

2^ 



6 = — 1 or 



(l + a) 2 
6(1 + a) + 1 



(3.6) 



(3.7) 



(3.8) 



3(1 + a) 

and both values lead to C = from eqn (3.4). Since eqn (3.3) is satisfied exactly, 
there are no approximate solutions with C ^ 0. 

The choice b = — 2/[3(l + a)] corresponds to the flat Friedmann model. In 
this case, eqns (2.17) and (2.18) are satisfied if 



x 



(a-l)/a 



3/5 2 (l + a) 2 ' 1 ° ° (l + 3a) 2 ' 

and one can choose x Q — S Q — 1 providing one scales the r and t coordinates such 
that 

P= ^ , 7 = ^4 (3.10) 
1 V3(l + a)' (l + 3a) 1 J 

x = z -2a/(l + a)^ g = z -2/[3(l+ a) ] (3 n) 



9(1 + a) 



(3.9) 



This gives 



and the metric becomes 

ds * = 01 dt 2 _ 7 -2 z -4/[3(l+a)] ^2 _ r 2(l+3 a )/[3(l+a)] ^[3(1+0,)] dQ 2 (3^ 
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One can put it in a more familiar form by making the coordinate transformation 

t = p t , f = /3 -2/[3(l+a)] r (l+3a)/[3(l + a)] > (3^3) 

which gives 

rfs 2 = d p_ £4/[3(l+a)] ^2 + f 2^2] (3^ 

This is just the flat Friedmann solution with p = a/i. We also have 

/i= _l_ ; y = ^±0^ z (l+Sa)/[S(l + a)]^ M= l J ,2(l+3a)/[3(l+a)]_ (3^5) 

The choice 6 = — 1 corresponds to the self-similar Kantowski-Sachs (KS) 
model. This is compatible with eqn (3.7) providing 

" 2 - ~ 0>/ °' (316) 

/,From eqns (2.17) and (2.18) we also require 

o2 -(l+a)/a _ ^ /o 1? \ 

-(l + a) 2_ 4a 2- ^ i7 J 

Eqn (3.17) shows that we cannot take x c = S Q — 1 in this case but both x Q and 
S are determined in terms of a and /3. The constant 7 is not constrained at 
all. If we take (3 and 7 to have the values given by eqn (3.10) for a < and % 
times those values for a > 0, so that we have the same r and t scaling as in the 
Friedmann solution, then eqns (3.16) and (3.17) give 

X0= \W\J ' 5 ° 2= (l + 3«)(l-a) (sRJ • (3 - 18) 

[Carr & Koutras (1993) do not incorporate the % factors for a > but this is 
a less sensible normalization since it allows the metric to be complex.] We now 
have 

S=S z-\ x = x z- 2a/{1+a) (3.19) 
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and the metric is 

ds 2 = [3 2 x 2 dt 2 - r^Ts^z^^dr 2 - S 2 t 2 dn 2 . (3.20) 

The t coordinate is spacelike and the r coordinate is timelike for a > because 
of the % factors in (3 and 7. For — 1/3 < a < 0, t and r have their usual 
interpretation but, from eqn (3.18), the circumferential coordinate is timelike 
since S Q is imaginary. One can put the metric in a more familiar form by making 
the coordinate transformation 

t = (3x t, f = 7- 1 ((3x ) 2a/{1+a) x 1/a S 2 r , (3.21) 

which gives 

ds > = d p _ t ^/(i +a)df .2 _ (Su/PxofPdQ 2 . (3.22) 
This corresponds to a p = a/i KS solution. We also have 

2a 2 

M= (a -l)(3 a+ l) " Afa - (3 - 22) 
V is negative for < a < 1 (corresponding to tachyonic solutions), while M is 
negative for— 1/3 < a < 1 (corresponding to negative mass solutions). Solutions 
with a < —1/3 are therefore more physical in that V and M are positive; although 
such solutions have negative pressure and violate the strong energy condition, the 
required conditions might arise naturally in the early Universe. Note that eqn 
(2.4) does not apply in this case because there is no well-defined origin; eqn (3.19) 
implies that R is independent of r, so everything is on a shell. Instead the value 
of m must be interpreted as the mass of the whole Universe at time t, as indicated 
by eqn (2.5). 

• q = 0, B + C = 0. In this case, one can show that the only consistent 
solution is 

a = b = 0, (3.24) 



19 



i.e. x and S are constant. [The condition q = permits another value of b but 
this leads to negative C, so the condition B + C = cannot be satisfied.] Eqn 
(3.24) implies that A is zero and hence eqn (3.3) is satisfied identically, so there 
are no approximate solutions with A^O. The condition B + C = also requires 

Si = ^xt a)/2a (3.25) 



for a / 0. This corresponds to the exact self-similar static solution, with the 
metric being given by 

ds 2 = (3 2 x 2 z 4a ^ 1+a Ut 2 - r^TS^dr 2 - r 2 S 2 Q dn 2 . (3.26) 

This can be put in an explicitly static form 

ds 2 = r^l^dt 2 - ^xl'oSfdr 2 - r 2 dQ 2 (3.27) 

by introducing the variables 

f = rS , t= fl±^V^ ~ 2a/(1+a) ^ (1_a)/(1+a) - (3-28) 
VI — a/ 



The other relevant functions are 

2a 



= x -(l+a)/a r -2^ y = x -(l-a)/2a z (l- a )/(l + a)^ j^f _ __J±__ (3,29) 



If (3 and 7 have the same values as in the Friedmann solution, eqns (2.17), (2.18) 
and (3.25) imply that x Q and S are given by 

9(1 + 300(1 + 6a + a 2 )] 2a/{1+3a) c \(1 + 6a + a 2 )^/ 2 {a + 1/3) 1+Qn 1/(1+3q) 



(18a) 3 / 2 



2a 

(3.30) 

so there is only one static solution for each equation of state. Note that eqn 
(2.20) implies that V R = in this expected. It should be stressed that 

there is no static solution in the dust case, essentially because one cannot put 
the x/a term in eqn (2.15) to zero. [If one puts a = in eqns (3.29) and (3.30), 
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one obtains M = and the spatial metric components all diverge; also eqn (3.25) 
cannot be satisfied.] The self-similar static solution was first discussed by Misner 
& Zapolsky (1964) and subsequently by Henriksen & Wesson (1978) and Carr & 
Yahil (1990). 

Note that there is an interesting connection between the static and KS so- 
lutions: if one interchanges the r and t coordinates in metric (3.26) and also 
changes the equation of state parameter to 

a ' = -TT25- (3 ' 31) 

one obtains the KS metric (3.20). For a static solution with a normal equation of 
state (1 > a > 0), a' must lie in the range —1/3 to 0, so some negative pressure 
KS solutions are related to positive pressure static ones. However, KS solutions 
with — 1 < a' < —1/3 correspond to |a| > 1 and so do not give physical static 
solutions. Note that a — a' only for a — or a — — 1. The mass of both the 
static and KS solutions tends to as a — > 0, although the solutions do not exist 
in the limit a = itself. 



• p = 0, A + C = 0. The condition p = implies 

1 (\ - a\ 1 (\ - a\ 

and eqn (2.19) then requires that V tend to the finite value 

v* = /rVM 1-0 ^- 2 . (3.33) 

The condition A + C = now implies 

V?(l-q) + 2a _ (l-a)(K 2 + a) 

(V? -!)(! + «)' 2a(l + a)(V?-l)' 1 J 



while eqn (3.5) yields 



9 =sa^- (3 ' 35) 
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This is only a consistent solution of eqn (3.3) for large z if V 2 > 1 and for small z 
if V 2 < 1. However, eqn (2.17) implies that the latter condition leads to negative 
values of M (and hence to unphysical solutions) unless V 2 < a. Since V 2 > 1, 
eqn (3.34) implies that both a and b are positive, so the density goes to zero and 
the scale factor goes to infinity in these solutions. It should be emphasized that, 
since B ^ from eqn (3.4), this case does not lead to an exact solution. Eqn 
(2.17) gives 

M ~ z -[K 2 (i-a)+i+3a]/(K 2 -i)(i+«) ( 3>36 ) 

and this necessarily tends to zero as z — > oo. On the other hand, eqn (2.18) 
implies 

M - - ~ 2 [K 2 (l-a)-«(l+3a)]/(^-l)a(l +Q ) [fe 2 (K 2 _ 1} _ 2& _ 1]. ( 3 .37) 

2 

If the exponent of z in this expression is positive, M — > as z — > oo only if the 
term in square brackets does and this requires 6 = 1/(K — 1). Eqn (3.34) then 
gives a quadratic equation for V*: 

(1 - a)V? - 2«(1 + a)K - a(l + 3a) = (3.38) 

with the real positive solution 



a(l + a) + i/afa 3 — a 2 + 3a + 1) 

K = — y — -. (3.39) 

1 — a 

[This assumes that a, V* and z are all positive; for a < one would need to take 
the negative square root in eqn (3.39).] Note that eqn (3.38) implies that the 
exponent of z in eqn (3.37) is indeed positive (as assumed). 

The value of V # given by eqn (3.39) exceeds 1 only for a > 1/5, so these 
solutions do not exist in the dust case. [For a < 1/5 there are solutions as 
z rightarrowO with 1 > V 2 > a but these have negative mass.] Eqns (3.33) and 
(3.39) impose a relationship between x and S , so these solutions are described 
by just one independent parameter. Requiring that the right-hand-side of eqn 
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(3.37) tends to — 1/2 as z — > oo merely determines the second order terms in the 
expansions for x and S. The metric has the asymptotic form 

ds 2 ~ z 2V * 2 W* 2 -Vdt 2 - z^-^dr 2 - rW^-Vdn 2 (3.40) 

and this can be reduced to the Minkowski form with a suitable change of coor- 
dinates. Note that although MS — > as z — > oo, it does so slower than z -1 , so 
that the mass itself (m = rMS) diverges. 

The forms of V(z) for the Friedmann, static and KS solutions are shown in 
Figure (la) for the general a case and in Figure (lb) for the a — 1/3 case. The 
asymptotically Minkowski solutions for a > 1/5 are not included since they are 
not exact (viz. Minkowski has no matter). Note also that the Minkowski solution, 
although static, is distinct from the exact self-similar static solution given by eqn 
(3.26). A rather peculiar feature of the similarity solutions, which arose in the 
context of the KS model, is that the mass can go negative. This may seem 
unphysical but - in the context of the Big Bang model - Miller (1976) has given 
a possible interpretation of this in terms of "lagging" cores. In the a = 1/3 case 
(and only this case), one can show that there is a curve in the V(z) diagrams 
where M = 0: 

V 3 = -y/3/2^ 2 (V 2 -l/9). (3.41) 

and this is also shown in Figure (lb). One sees that the curve has asymptotes 
at V — ±1/3. The upper part (with V > 1/3) is relevant for asymptotically 
Friedmann solutions, while the lower part (with V < —1/3) is relevant for asymp- 
totically KS solutions. M is negative in between the two parts and this region 
includes KS itself (as expected). Note that eqn (3.41) is not a sufficient condition 
for M = 0; in fact, it implies that M has two possible values, only one of which 
is zero. 
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3.2 Logarithmic Power-Law Similarity Solutions 



By analogy with eqn (3.1) we now look for solutions in which z tends to some 
finite value z* and in which 



x = x L a , S = S L b , L = ln{z/z,) 



(3.42) 



for constants xq and Sq. Clearly z = z* corresponds to an infinite distance from 
the origin for b < and zero distance for b > 0. Eqn (2.15) requires 



6-1 + {l + a) 



2b 



1 + a c 

and the leading term is zero only for 



b= 3 + 



+ 



1 + a 
3a 



3 + a 
1 + a 



b- a - 

a 



L = 



a. 



(3.43) 



(3.44) 



It turns out that the last term is never zero, so these are only asymptotic and 
not exact solutions. There are now two possible situations, according to whether 
V tends to infinity or some constant value V*. 



V — > V* as z — > z*. In this case, eqn (2.16) can be written in the form 



! _ 2V?ab x 2a 2 q(l + q)/gVg 2(1 . tt)/(1+a) ( , r w«-i)/« 

~(V?-a) L + (l + a)(K 2 -«) K 2 "« * ( ° j 

(3.45) 

The only consistent solution to this equation has the last term tending to zero 
and this then implies 

/ 2V 2 a \ 

X b. (3.46) 



K 2 - a 



However, eqn (2.19) also implies 

2a 



so we require V 2 = 1. Eqns (3.44) and (3.47) determine a and b and lead to 

S = S £ (1_Q)/(1 ~ 5a) , x = a;o J L 2Q/(1 ~ 5Q) . (3.48) 
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Thus the scale factor diverges and the density goes to zero providing a > 1/5. 
However, there are no such solutions for a < 1/5, so they do not appear in the 
dust case. The condition V* — 1 gives a relationship between the constants x , 
So and z*, so these solutions are described by two independent parameters. Eqn 
(2.17) implies 

M ~ x,( 1-a )/( 5a-1 ) ) (3.49) 

so this is zero at z = z*. However, eqns (3.48) and (3.49) imply that the product 
MS tends to a constant, so the mass itself is not zero. Eqn (2.18) can be written 

as 

M = \ + \?x-*l°& | (| ) 2 (V 2 - 1) - ^ - 1 1 . (3.50) 
Since x~ 2 ^ a S 6 — > oo, this can go to zero only if 

f(V 2 -l)-2 (3.51) 



and eqn (2.19) then implies 

V 1 -5a 



V 1-a 

However, since the last term in eqn (3.50) scales as 



L (a-l)/(5a-l) 



(3.52) 



(3.53) 



we also need the term in square brackets to go to zero as l/ 1-a )/( 5a-1 ). This 
determines the second order terms in the expressions for x and S but does not 
impose any further relationship between xo, So and z*. Note that the metric can 
be written as 

ds 2 ~ L 4a ^-^[dt 2 - dr 2 - r'L^-^^dQ 2 } (3.54) 

This resembles the open Friedmann model and can be transformed to Minkowski 
form with a new choice of time-slicing (as in the Milne model). These solutions 
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are therefore asymptotically flat or, more precisely, asymptotically Schwarzschild 
(since the mass tends to a non-zero constant). They are discussed further by Carr 
et al. (1998). 



• V — > oo as z — > z*. In this case, eqn (2.16) can be written in the form 

aL- 1 = lahL- 1 - a(l + a)(3 2 z; 2 ^^ 1+ ^ (x £) a(Q ~ 1)/a . (3.55) 

It is easy to show that the only consistent solution has a = a/(l — a), so that all 
the terms scale as L -1 , and eqn (3.44) then implies b = 2/[3(l — a)], so we have 

S = S L 2 / 3(1 - a \ x = x L a ^ 1 - a \ (3.56) 

Eqn (2.17) then gives 

M~L- 2/3(1 - a) , (3.57) 

so M — > oo and MS tends to a constant at z — z*. Note that eqn (3.55) yields the 
same relation between A, B and z* as eqns (2.17) and (2.18), so these solutions 
are described by two independent parameters. Eqns (2.13) and (2.14) imply that 
the metric tends to 

ds 2 ~ L 2a /^dt 2 - L-W^dr 2 - r 2 L 4 / 3 ^dn 2 : (3.58) 

corresponding to a Schwarzschild-type singularity of infinite density. 

4 Self-Similar Solutions with Dust 

We first discuss the dust similarity solutions because these exhibit many of the 
features of the similarity solutions with pressure, even though the equations are 
considerably simplified. In particular, if a = 0, there is no sonic point and both 
the mass and energy within comoving radius r are conserved, so E and m/r = MS 
are constant. If we put m/r = k, then eqn (2.4) implies 

o<9i? dm 
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and this can be combined with eqns (2.12) and (2.14) to give 



e x = 7 - 1 (47r/ir 2 5 



2 o2\-l 



k 7" 



dr 



(4.2) 



On the other hand, the definition of T in eqn (2.6) implies e x = T 1 (dR/dr) and 
so the constant k is just r/7. The mass function is therefore 



M 



V1 + 2E 



(4.3) 



7 5 7^ ' 

where we have taken the positive square root to ensure that the mass if positive. 
(We discuss the negative mass case later.) Since eqn (2.13) implies e u = j3, eqn 
(2.6) can now be rewritten as 

2 VT+2E 



E = 



2(3 2 



dz 



IS 



(4.4) 



This equation is equivalent to eqns (2.17) and (2.18) with a — 0. Indeed, from 
eqn (2.12), providing one makes the substitution 

1 x dln(nr 2 ) 



x 



1, 



(4.5) 



a x dlnz 
eqns (2.15) to (2.18) with a = are all formally satisfied. 

In this case, rather than imposing eqn (3.10), it is convenient to scale the r 
and t coordinates so that f3 — 7 = 1. This is equivalent to making the rescaling 
r — > 7r and t — > t/(3. Eqn (4.4) then implies 



dS_ 

dz 

and this can be integrated to give 



2E + 2TIS 



(4.6) 



VES 2 +TS 2T ^ n U-l [ES 

V2E (2£)3/ 2 611111 V T 



D±- 

z 



yJ 5 3/2 



VES' 2 +rs 



V2E 



± 



2F ■ -1 

■ sin 



(-2£)3/ 2 



(E > 0) 
(E = 0) 
p (-1/2<S<0) 



(4.7) 
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where D is an integration constant, the plus and minus signs apply for dS/dz 
negative and positive, respectively, and sin -1 is always taken to be positive. If 
we took the negative square root in eqn (4.3), corresponding to M and Y being 
negative, there would be another solution for E > given by 



1 JES 2 -\T\S 2|r| 1 [ES 

This solution is unphysical, since the mass is negative, but it is of interest for 
later comparison with the solutions with pressure. Eqns (2.10) and (2.11) give 
the velocity functions as 



Sz ± J2E + 2Y/S 
V = V- '- (4.9) 



and 



J2E + 2T/S 
Vr = ±^ fr ^- 1 (4.10) 



while eqns (2.10) and (4.2) imply that the density is given by 



1 V1 + 2E 



47T/Xt 2 = — — = (4.11) 

zS 2 V zS 2(Sz ± pE + 2Y/S) 

where the plus and minus signs again correspond to dS/dz being negative and 
positive, respectively. Note that V and \i are negative (corresponding to tachyonic 
models) for the solution given by eqn (4.8). 

Eqn (4.7) implies that there is a 2-parameter family of similarity solutions (as 
in the general a case). The flat Friedmann solution has D = E = 0, so we need 
to determine the significance of the solutions in which D and E are non-zero. In 
each case, we will show the form of S, V and fit 2 as functions of z. In obtaining 
the full family of solutions, it is crucial that we allow z to be either positive 
or negative. Our analysis will also cover the unphysical solutions with negative 
mass. 
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4.1 D = solutions 



Solutions with D = are asymptotically Friedmann as \z\ — > oo and are specified 
entirely by the energy parameter E 1 . The forms of S^z) and V(z) in these solutions 
are shown in Figures (2a) and (2b). The solutions with z > correspond to 
initially expanding Big Bang models and have V > everywhere: they start 
from an initial Big Bang singularity (S — 0) at t — (z — oo) and then either 
expand indefinitely (S — > oo) as t — > oo (z — > 0) for E > or recollapse to a 
black hole singularity (S = 0) at 

for _E < 0. In the first case, V decreases monotonically from oo to 0. In the 
second case, it reaches a minimum before rising to oo at z s ; the minimum will 
exceed 1 (in which case the whole Universe is inside the black hole) if E is less 
than some critical negative value E* and it will be less than 1 (in which case there 
is a black hole event horizon and a cosmological particle horizon) if E exceeds 
E*. Such solutions were studied in detail by Carr & Hawking (1974) and Carr 
& Yahil (1990). The solutions with z < are the time-reverse of these and have 
V < everywhere: as t increases from — oo to (i.e. as z decreases from to 
— oo), the E > models collapse from an infinitely dispersed state (S = oo) to a 
Big Crunch singularity (S — 0); the E < models also collapse to a Big Crunch 
singularity but they emerge from a white hole and are never infinitely dispersed. 

Both S and V have the same z-dependence as in the E = Friedmann solution 
as \z\ — > oo: 

S = [9Vl + 2E/2] 1/3 \z\- 2/3 , V = [6(1 + 2E)]- 1/3 z 1/3 . (4.13) 

However, the E ^ solutions deviate from the E — solution at small values of 
\z\. The E < solutions never reach z = at all, while the E > ones have 

S =(2E) 1/2 \z\-\ V = -(l + 2E) 1/2 E- 1 zln[(2E) 3/2 (l + 2E) 1/2 \z\}. (4.14) 
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as \z\ — > 0. The first relation implies that the circumference function R(r, t) = Sr 
is non-zero in limit r — > unless = since 



This means that the "coordinate" origin (r = 0) is an expanding 2-sphere. This 
has a natural physical interpretation since the forms of S and V are similar to 
those in the KS solution [cf. eqns (3.19) and (3.23)], in which all the matter is 
localized on a shell. [However, there is no exact KS solution in the dust case.] 
To obtain a complete solution, one must match the self-similar solution onto a 
(non-self-similar) void inside R(t, 0). In the E < case, the physical origin is the 
black hole singularity (since R = rS = there), so only in the E = case can 
one identify z = with the physical origin. 

The form of the density function fit 2 is shown in Figure (2c). At a given time 
this specifies the density profile fx(r) and it illustrates that a non-zero value of 
E necessarily introduces inhomogeneity into the model. Solutions with E > 
are everywhere underdense relative to the Friedmann model, whereas those with 
E < are everywhere overdense. For the latter, the density diverges at the 
singularity and eqns (4.3) and (4.12) imply that the mass associated with this 
singularity is 



It therefore starts off zero when the singularity first forms at t — but then grows 
as t. 

The form of the mass function M(z) in the D = solutions is not shown 
explicitly but can be immediately deduced from the expression for S since eqn 
(4.3) implies M ~ S -1 . In the E > case, there is always a single point where 
M — 1/2 and this corresponds to the cosmological apparent horizon. In the 
E < case, eqn (4.7) implies that S has a maximum oiT/\E\ and so eqn (4.3) 
shows that M has a minimum of \E\. Since this is less than 1/2, there are always 




(4.15) 



m s = (MSz) s t = {-2E) 3/2 t/{2n). 



(4.16) 
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two points where M — 1/2, one corresponding to the black hole apparent horizon 
and the other to the cosmological apparent horizon. Note that a black hole's 
apparent horizon always lies within or coincides with its event horizon, which is 
why the first can exist without the second. 

Finally, we consider the E > negative-mass solutions given by eqn (4.8). 
Their form is indicated by the broken curves in Figures (2). S, V and /i have the 
same form as in the positive mass solutions for small values of \z\ except that V 
and \x reverse their signs. However, the solutions are very different at large values 
of \z\ since eqn (4.6) shows that S must always exceed |r|/i£. Indeed it tends to 
this value asymptotically, so we have 

S « Vl + 2E/E, Vm-z/E, /2r 2 ^ (4.17) 

as \z\ — > oo. The form of this solution is closely related to that of the a « 1 
static solution [cf. eqn (3.29)], although there is no static solution in the a = 
case itself. 

4.2 E = solutions 

We now put E — and consider the effect of introducing a non-zero value for 
the constant D. [In this case, eqn (4.6) does not permit T < 0, so there are no 
negative-mass solutions.] The form of S(z) for the D > solutions is shown in 
Figure (3a). There are two types of solutions in this case, one expanding and the 
other collapsing. For the expanding solutions (solid lines), S = at z = — l/D 
and so the Big Bang occurs before t — (i.e. it is "advanced"). As t increases to 
(i.e. as z decreases to — oo), S tends to the finite value 

Soo(D) = (3D/V2) 2 / 3 . (4.18) 

As t further increases from to +oo (i.e. as z jumps to +oo and then decreases 
to 0), S increases monotonically to oo. For the contracting solutions (broken 
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lines), S starts infinite at t — — oo (z = 0) and then decreases to S^D) as t 
increases to (z — > — oo). As t further increases (i.e. as z jumps to +00 and then 
decreases), S continues to decrease until it reaches at the Big Crunch singularity 
at z — 1/D. Both types of solutions are characterized by the fact that they have 
just one singularity. 

The form of V(z) in the D > solutions is shown in Figure (3b). For the 
expanding solutions (solid lines), it starts off at —00 at the Big Bang (z = —1/D), 
reaches a negative maximum (which will be more than -1 if D exceeds some critical 
value D*) and then, from eqn (4.9), tends to 



as z — > —00. When z jumps +00, V becomes positive but eqn (4.18) still applies. 



any solution with finite S at infinity must be "nearly" static in the sense that 
dS/dz tends to zero. However, the solutions are not asymptotic to an exact static 
solution (indeed this does not exist in the a = case) because eqn (4.10) implies 
that Vr tends to a non-zero value: 



This is also reflected by the fact that, in terms of the variable £ = 1/z, dS/d£ 7^ 
at £ = 0. We therefore term these solutions asymptotically "quasi-static". As 
z decreases from +00 to 0, V decreases monotonically to 0. For the contracting 
models (broken lines), V starts from zero at z = and monotonically decreases as 
z decreases to —00, being given by eqn (4.19) asymptotically. When z jumps to 
+00, V jumps to +00 and then decreases to a minimum before rising to infinity 
at the Big Crunch singularity. The minimum will be less than 1 if D exceeds the 
value -D*; in this case, one necessarily has a naked singularity. 

The form of the density function fit 2 in the D > solutions is interesting. 



V = S OQ (D)z=(3D/V2) 2 / 3 



z 



(4.19) 



For large \z\ : V ~ z rather than z 1 ^ 3 (as applies in the D = case) because 




(4.20) 
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From eqn (4.11) the density parameter has the form 

where the + and - signs apply for negative and positive values of dS/dz, re- 
spectively. (The factor Qirt 2 corresponds to the density in a flat Friedmann dust 
universe.) For a given fluid element, this just describes how the density evolves 
as a function of time and it has the expected form. However, at a given time 
it also prescribes the density profile and one sees immediately that a non-zero 
value of D (like a non-zero value of E) introduces an inhomogeneity. This inho- 
mogeneity has a particularly interesting form. In the z < regime, the profile 
for the collapsing solutions is homogeneous for \z\ << 1/D but roughly static 
and isothermal (fi ~ r~ 2 ) for \z\ » 1/D. In the z > regime, the profile 
corresponds to a density singularity at the centre of a nearly static isothermal 
sphere. These features are illustrated by the broken lines in Figure 3(c) and have 
an obvious physical interpretation. It is interesting that the isothermal model 
features prominentl! y in both regimes despite the fact that there is no exact 
static solution in the dust case. The mass of the singularity in these solutions is 

m s = (MSz)s t = t/D (4.22) 

from eqn (4.3). As in the asymptotically Friedmann case, it therefore starts off 
zero at t — but then grows as t. As pointed out by Ori & Piran (1990), such 
solutions have an analogue in Newtonian theory (Larson 1969, Penston 1969). 

Finally, we consider the D < solutions. The form of S(z) in this case 
is shown by the dotted curves in Figure (3a). Such solutions are confined to 
\z\ < 1/D, with S either decreasing monotonically for z < (i.e. as t increases 
from — oo) or increasing monotonically for z > (i.e. as t increases to +oo). 
However, these solutions break down before S reaches 0. This is because, as 
indicated by the dotted curve in Figure (3b), \V\ increases to some maximum 
value and then falls to zero at \z\ = 1/(3D). From eqn (4.20) this means that 
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the density diverges there, as indicated by the dotted curve in Figure (3c). This 
divergence is associated with the formation of a shell-crossing singularity since the 
model resembles the KS solution at this point. For \z\ > —1/(3D), the density 
and velocity function become negative but this is presumably unphysical. 

4.3 D^O^^O solutions 

The forms of S(z) and V(z) for the (D > 0, E ^ 0) solutions are indicated by the 
solid and broken lines in Figures (4a) and (4b). The figures assume that D is 
fixed but allow E to vary. The (D > 0, E > 0) solutions are qualitatively similar 
to the (D > 0, E — 0) ones (i.e. one has monotonically expanding or collapsing 
models). In particular, the form of the solution near z = is still given by eqns 
(4.14), so the value of D is unimportant here. As before, z = is no longer the 
origin, so one has to attach the solution to a non-self-similar central region. As 
\z\ — > oo, the solutions with E > have 

S ~ Soo(D, E), V « S 00 (D,E)(l + 2E)- 1 '*z (4.23) 

where S^D, E) is the solution of eqn (4.7) with z = oo. As in the (D > 0, E — 0) 
case, these solutions are only quasi-static and not exactly static at large \z\ since 
eqn (4.10) implies 

V^(oo) = T- 1 y /2E + 2T/S 00 (4.24) 

and this is always positive. Also as in the (D > 0, E — 0) case, |V| reaches a 
minimum before diverging at the Big Bang or Big Crunch singularity and this 
minimum will fall below 1, corresponding to a naked singularity, providing E 
exceeds some critical negative value E*(D). The dependence of E* on D can only 
be determined numerically but it clear that this corresponds to lying outside some 
boundary in the (E, D) diagram. 

The (D > 0, E < 0) solutions are qualitatively different from the (D > 0, E — 
0) solutions in that there are no models which collapse from or expand to infinity. 
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This is clear from eqn (4.6) which implies that S never exceeds One might 

expect that there could be an exact static solution with S = T/(—E) but this is 
incompatible with eqn (2.15) if a = 0. The models which expand from the Big 
Bang singularity at z — — 1/D recollapse to a black hole singularity at 

=-1 -l 



zs = 



D 2irVT+2E' 



(4.25) 



-2£) 3 /2 

(necessarily less than 1/D), whereas the models which recollapse to the Big 
Crunch singularity at z = 1/D emerge from a white hole singularity at —z s . 
Again there will be a naked singularity if the minimum value of |V| is less than 
1. 

The form of fit 2 in these solutions is shown in Figure (4c). It is very similar to 
the form indicated in Figure (3c) except that one has added curves corresponding 
to the E < solutions. The E > solutions are everywhere underdense relative 
to the E = solutions, whereas the E < solutions are everywhere overdense. 
The form of M(z) can be deduced immediately from Figure (3a) and eqn (4.3). 
As \z\ — > oo, M — > Vl + 2E/S<x>(D, E) and this exceeds 1/2 for sufficiently low 
values of D. In this case there is no black hole or cosmo logical apparent horizon, 
which is different from the D = case. 

The D < solutions have the same form as in the E = case, except that the 
E < solutions emerge from a white hole rather than collapsing from infinity. 
As in the E = case, such models are probably physically unrealistic since the 
density diverges due to shell-crossing. They are therefore not shown explicitly. 
The form of the (unphysical) negative-mass solutions, which only exist for E > 0, 
is indicated by the dotted curve in the Figures (4) and is similar to the D = 
case shown in Figure (2). 
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5 Self-Similar Solutions with Pressure 



One might expect the a^O solutions to share some of the qualitative features 
as the a = ones in the supersonic regime. The arguments for this are partly 
physical (viz. pressure effects should be unimportant on sufficiently large scales) 
and partly mathematical (viz. the dust equations can be obtained from the 
general a equations by taking the limit a — > 0). The behaviour of the self-similar 
dust solutions therefore suggests that there should exist at least two classes of 
similarity solutions with pressure at large \z\: a 1-parameter family which are 
asymptotically Friedmann and a 2-parameter family which are asymptotically 
quasi-static. In this section we will show that there are indeed solutions of these 
kinds. However, we saw in Section 3 that new possible behaviours arise at large \z\ 
when there is pressure. In particular, there is an exact static solution and an exact 
Kantowski-Sachs solution, so one might expect there to be families of solutions 
asymptotic to these. We will demonstrate that this is indeed the case and that 
each of the families is described by one parameter. We saw in Section 3 that 
there are also asymptotically (but not exact) Minkowski solutions for a > 1/5. 
We discuss these solutions in more detail elsewhere (Carr et al. 1998) but do not 
consider them further here. 

The inclusion of pressure also introduces qualitatively new features in the 
subsonic regime, so there are important differences from the dust solutions at 
small values of \z\. In particular, the presence of a sonic point at |V| = \fa 
allows solutions to be discontinous there, so one might anticipate a wide variety 
of subsonic behaviours. However, the requirement that the solution be regular at 
the sonic point severely restricts the possible subsonic behaviours, with the con- 
sequence that the only possible solutions with pressure at small \z\ are the exact 
static model, a 1-parameter family of asymptotically Friedmann models and a 1- 
parameter family of asymptotically Kantowski-Sachs models. Indeed the value of 
z at the sonic point almost uniqely (modulo the number of oscillations) identifies 
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the subsonic solution, which means that one can effectively extend a solution in 
the supersonic regime into the subsonic regime in only one way. Although other 
extensions exist, they will become unphysical at some point (either because the 
mass becomes negative or because they encounter another sonic point at which 
they become irregular). 

We will start by discussing the characteristics of the asymptotically Friedmann 
and asymptotically KS solutions. In this context we will be mainly reviewing the 
work of Carr & Yahil (1990) and Carr & Koutras (1993) but we will extend these 
earlier studies somewhat, so as to make the connection with the dust solutions 
more explicit. We will then discuss the asymptotically quasi-static solutions. 
The discussion here will be entirely original, although some examples of this type 
of solution have been considerd before. The form of the function V(z) for the 
various solutions is shown in Figures 6-8. The a = 1/3, a > 1/3 and a < 1/3 
cases are shown separately, since they have somewhat different characteristics. 
The a = 1/3 case is special since the equations simplify, so this is the one which 
has been most studied numerically. To avoid an unnecessary proliferation of 
detail, we only present the z > curves. The z < curves are trivially derived 
from these by reflecting in the origin (cf. Figures 2-4): one merely replaces z by 
\z\ and V by —V. However, it is important to stress that, as in the dust case, 
the complete solutions may involve both the z < and z > regimes. 

5.1 Asymptotically Friedmann Solutions 

Carr & Yahil (1990) consider solutions which are either exactly or asymptotically 
Friedmann for large and small values of z. They therefore introduce functions 
A(z) and B(z) defined by 

x = z- 2a ^e A , S = z- 2 ^ 1+a ^e B . (5.1) 
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The Friedmann solution itself (A = B = 0) passes through the sonic line Q at 

■ -, 3(l+a)/(l+3a) 

1 + 3a 

For comparison the static solution passes through Q at 



(5.2) 



zs 



/2o;)3/23(5a-l)/2(a-l) "1 (l+a)/(l+3a) 



(5.3) 



(1 + 3a)(l + 6a + a 2 )_ 

and one can show that this is always less than zp. In fact, the only physical 
subsonic solution which passes through zp is the Friedmann solution itself. All 
the other solutions are unphysical because, as z decreases, V either reaches a 
minimum and then hits the sonic surface again but off the line Q or the mass 
within r goes negative (Bicknell & Henriksen 1978a). Since the only solution 
which is exactly Friedmann outside the sonic point is the Friedmann solution 
everywhere, we confine attention to solutions which are asymptotically Friedmann 
at large z. 

The ordinary differential equations for x and S now become ordinary differ- 
ential equations for A and B. If we linearize the equations in A, A and B to find 
the 1st order solution as z — > oo, eqns (2.15) and (2.16) yield 



B = ( -L ] 4 



B 



1 



A 



l + 3a 
1 + a 

(a -I) 



B, 



-A 



(5.4) 



(5.5) 



,2a/ 3a(l + a) 

Differentiating the second equation and eliminating B and B then leads to the 
following differential equation for A: 



A + 



This has two solutions, 



(9a - 1)^ 2(1 + 3a)(a 



3(1 + a) 



3(1 + 



0. 



(5.6) 



A 2-2(l+3a)/3(l+a) Qr A ^ z (l-a)/(l+a) 



(5.7) 
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but the second one can be rejected since the exponent is positive for a < 1 (so 
that A diverges as z — > oo). The first solution then gives 

A = _ a ( 1 + 3 ^) ^-2(l + 3a)/3(l +a ) ? B = Boo - ^- 2 d+3°)/3(l+«) ( 5 . 8 ) 

(1 + a) 

where Boo and /c are integration constants. Note that A — > as z — > oo because 
eqns (2.12) and (5.1) show that A has physical significance (viz. the density 
perturbation) but there is no physical restriction on (the asymptotic value of 
B). The constants in eqn (5.8) are related since eqns (2.17) and (2.18) imply 

3(l + a)(e-2*--e**~) 

2(l + 3a)(5 + 3a) ' 1 ' 

Thus there is a 1-parameter family of asymptotically Friedmann solutions. [De- 
spite the presence of the parameter B^, these solutions are really asymptotic to 
the exact Friedmann model, since one could formally gauge B^ to zero at infinity 
(but not finite z) by taking a different spatial hyper surf ace.] /,From eqn (2.20), 
the energy function is 

E = E OD + (z-^ 1+3a)/ ^ 1+a) ), = -(e 6B ™ - 1) (5.10) 

2 

Eqn (2.12) shows that the solutions are overdense or underdense relative to the 
Friedmann solution according to whether A < or A > 0, respectively. From eqns 
(5.9) and (5.10), this corresponds to (k > 0,5^ < 0, E^ < 0) or (k < 0,5^ > 
0, Eoo > 0), respectively. In the Friedmann case itself, k = B^ = E^ = 0. 

If B^ is sufficiently negative, V reaches a minimum value Knin above ^fo. 
as z decreases and then rises again to infinity. Such solutions are supersonic 
everywhere and contain black holes which grow as fast as the Universe. There is 
an event horizon and particle horizon providing V m [ n < 1 and this will apply if 
B^ is more than some critical negative value B^; otherwise the whole Universe 
is inside the black hole. However, since eqn (2.17) implies that M = at S = 
and eqn (2.18) then implies M < 1/2, there is still an apparent horizon for 
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-Boo < B^.; this generalizes the result found in the dust case. As increases (i.e. 
as the asymptotic overdensity decreases), the values of V m i n and z m { n decrease. 
Eventually it reaches another critical negative value B^ %t at which V min = y/a 
and, for > B^ lt , the solutions must reach the sonic surface. /,From eqn 
(5.10), these critical negative values for B^ also correspond to critical negative 
values for the asymptotic energy, E* and E crit . As B^ continues to increase, the 
value of z at which the solution goes transonic (z s ) increases, passing through the 
value indicated by eqn (5.2) when B^ = and tending to infinity as B^ goes to 
infinity (corresponding to increasingly underdense solutions). All the solutions 
with B^ > B^ lt reach the sonic surface but only the ones which cut it on the 
line Q are regular. This only applies if z s lies within the ranges Z\ to z 2 or above 
z 3 [as indicated in Figure (1)]. A diverges at the sonic point for values of B^ 
corresponding to z 2 < z s < z 3 . 

We next consider the condition that the solution be well-behaved at z — 0, in 
the sense that the density and velocity are finite. This requires that A and B be 
finite, which implies 

V(0) = 0, A(0) = B(0). (5.11) 
We also require M(0) = 0, and eqn (2.18) then implies 

A(0) = 3aS(0), (5.12) 

which shows that there is a 1-parameter family of solutions which are well-behaved 
at the origin. We will take this parameter to be A Q = A(0). This is a measure of 
the overdensity at the origin since, from eqns (2.12) and (5.1), 

>f(0)~ 



A 



a 
1 + a 



log 



Mo) 



(5.13) 



where /j,f is the density in the Friedmann solution. Thus A > and A < 
solutions are underdense and overdense, respectively. Eqn (2.6) implies that 
E — > as z — > 0, so the energy function goes to zero as S — > oo (as expected). 
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For some range of values of A the subsonic solutions must hit the sonic surface 
in (x, S, S) space, since the solution with A = does. As before, the regular 
solutions must hit it on the line Q. However, the behaviour of the subsonic parts 
of the solutions is more complicated than that of the supersonic solutions. As 
the parameter A decreases from positive values to some critical negative value 
Aq™*, z s decreases continuously to z\. In this parameter range the solutions with 
z s > z 3 are regular at the sonic point, while those with z\ < z s < z 3 are all 
irregular. [In the a — 1/3 case, for which z 3 = Zp, one therefore has a continuous 
family of underdense solutions but no overdense solutions within this range.] As 
Aq decreases below Aq"*, the V(z) curves develop an inflexion and z s increases 
again to the value z 2 . The subsonic solutions thus cross over each other in V(z) 
space. Although one does not reach every value of z s between z\ and z 2 , there 
is a band of solutions within z\ < z < z 2 which are regular at the sonic point. 
This corresponds to the first band of overdense solutions and is associated with 
just a small range of A values. As A decreases further, z s moves back and forth 
between the values Z\ and z 2 and the V(z) curves exhibit an increasing number of 
oscillations. One can group these solutions into families according to the number 
of oscillations they exhibit. Each family will contain a narrow band of solutions 
(corresponding to bands of A Q values) which are regular at the sonic point. The 
existence of these bands was first pointed out by Ori & Piran (1990). Indeed this 
band structure arises even in the Newtonian situation, as shown by Whitworth 
& Summers (1985). Note that all the overdense solutions are nearly static close 
to the sonic point (in fact, z 2 = z s if a = 1/3), although they deviate from the 
static solution as one goes towards the origin. 

Provided there are points on Q which are intersected by solutions which are 
asymptotically Friedmann at both large and small z, one can construct a solution 
with a sound-wave which represents a density perturbation growing at the same 
rate as the Universe. One would expect this to be possible providing there is a 1- 
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parameter family of solutions at the sonic point, i.e. providing z s lies in the range 
of values between z\ and Z2 and above Z3. For in this case, for each point on Q, 
one would expect at least one supersonic solution to be asymptotically Friedmann 
and at least one subsonic solution to be regular at the origin. Furthermore, one 
would expect the value of A (corresponding to the density gradient or velocity 
gradient) to be continuous at the sonic point in such solutions since only one 
value corresponds to the 1-parameter family. The numerical calculations of Carr 
& Yahil (1990) show that transonic solutions do indeed exist and have the features 
aniticipated. In particular, the solutions all have a continuous velocity gradient at 
the sonic point, although they do not span the entire range of values Z\ < z s < Zi- 
Note that for each a there is also one asymptotically Friedmann solution which 
can be attached to the exact static solution inside the sonic point. 

The form of the asymptotically Friedmann solutions for a = 1/3, the case most 
likely to pertain in the early Universe, is shown in Figure (5). The qualitative 
form of the solutions for more more general values of a is indicated in Figures 
6-8. In deriving the form of the velocity function, we use eqns (2.19), (5.1), (5.8) 
and (5.12) to express the asymptotic behaviour as 

V « e {l - Za)A ° /a V F {\z\ « 1), V « e- 2B °°V F {\z\ » 1) (5.14) 

where V F is the exact Friedmann velocity. For a > 1/3, V is more (less) than V F 
for the overdense (underdense) solutions at all values of z. However, for a < 1/3, 
V starts below (above) V F at small z and ends up above (below) it at large z for 
the overdense (underdense) solutions. For a = 1/3, the dependence of V upon 
A Q at small \z\ only appears at second order. 

5.2 Asymptotically Kantowski- Sachs Solutions 

If we wish to consider solutions which are asymptotic to the self-similar KS model, 
then eqn (3.19) suggests that we introduce functions A(z) and B(z) defined by 

x = x Q z~ 2a/(1+a) e A , S = S z- 1 e B . (5.15) 
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As in the asymptotically Friedmann case, we linearize the equations in A, A and 
B to find the 1st order solution as |V| — > oo. Eqns (2.15) and (2.16) then yield 

If we differentiate the second equation and then substitute for B using the first, 
one obtains the following differential equation for A: 



[ a + lJ 1 + a) 2 



This has two solutions: A oc z Pl and A oc z P2 where 

-l + a± v /(l-a)(24a 2 + 7a + l) 

V ^ = 2(TT70 • (5 ' 19) 

For a > 0, the KS solution has |V| — > oo as z — > oo, so we must choose the 
positive root pi. The general solution then has the form 

1 n-a\ 1 



^4 — A^z P1 , B — A^ 



z- pi (5.20) 



2a VI + a J pi 

where A^ is an integration constant. For —1 < a < —1/3, the KS solution has 
\V\ — > oo as z — > 0, so we must choose a negative root. Only p2 is negative for 
this range of a, so this gives a solution like eqn (5.20) but with p2 replacing p\. 
In both cases there is thus a 1-parameter family of asymptotically KS solutions. 
For —1/3 < a < 0, one again has |V| — > oo as z — > oo but both p\ and P2 are 
negative, so there is no solution as z — > oo. [Since eqn (3.31) shows that KS 
solutions with — 1/3 < a < also correspond to static solutions with < a < 1 if 
r and t are interchanged (so that z goes to 1/z), this is related to the fact - shown 
later - that there are no asymptotically static solutions as z — > for < a < 1.] 
Eqn (2.20) implies that E — > — 1/2 as 2; — > 00, so the dust analysis suggests that 
these are in some sense extreme black hole solutions. 
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At small values of |V|, A = B = and eqns (2.17) and (2.18) imply that the 
values of A and B at z = are related by 

e 2B ° = l -[M KS e- A ^/« - {M KS ~ \)e~' A r\ (5-21) 

where M KS is defined by eqn (3.23). Again there is a 1-parameter family of 
solutions and A is a measure of the underdensity or overdensity at the origin 
relative to the exact KS solution. Note that there are only isolated solutions at 
a sonic point for < a < 1 and this means that any asymptotically KS solution 
which hits the sonic surface is unlikely to be regular there. For although there is 
a 1-parameter family of solutions from each point of the line Q, corresponding to 
the freedom in the value of A, only two of these are regular and they are both 
isolated. It is therefore unlikely that either of them will be asymptotically KS 
since the general solution has two parameters, whereas the asymptotically KS 
solutions have only one. 

Carr & Koutras (1993) have integrated the equations in the a = 1/3 case, 
although the physical significance of these solutions is unclear. Figures (6) show 
the form of V(z) in this case. Let us first consider the supersonic solutions. The 
underdense solutions have A^ positive. As z decreases, they all cross V = — 1 at 
some point to the left of the exact KS solution. However, they do not hit the sonic 
point but reach a maximum between V = —l/y/3 and V = — 1 as z decreases. 
They then hit the V = — 1 surface again (all with the same value of z), with M and 
/j, tending to and the scale factor S diverging. [This behaviour is analogous to 
that which arises for the solutions which are asymptotically Minkowski at finite 
z] The overdense supersonic solutions have A^ negative and, as z decreases, 
they all hit the sonic line to the right of the exact KS solution. As A^ decreases, 
the point at which they hit the sonic line moves to infinity. All the supersonic 
solutions have M < everywhere and so never cross the M = curve. Let us 
now consider the subsonic solutions. The overdense ones have A negative. None 
of the curves hit the sonic surface since they reach a minimum as z decreases and 
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then asymptotically approach V = 0. The interesting feature of these solutions 
is that the function M, which is negative at the origin, goes through [i.e. the 
solution crosses the curve M = in V(z) space] and eventually becomes positive 
as z increases. M and \i tend to infinity as z — > oo. The underdense solutions 
have A positive and hit the sonic line to the left of KS. 

5.3 Asymptotically Static and Quasi-Static Solutions 

If we wish to consider solutions which are asymptotically static, we introduce 
functions A(z) and B(z) defined by 

x = x Q e A , S = S Q e B (5.22) 

where x Q and S are given by eqn (3.30). Eqns (2.15) and (2.16) then become 

B + 3B 2 - i + (?±*\ B - ( l -±±\ AB = (5.23) 



a \a + 1/ V a 

V 2(b-A.) = -A + ( [e -4B+A(l-«)/a _ !]_ (524) 



la) 2 VI + a, 
To find the first order solution as V — > oo (i.e. as z — > oo), we linearize these 
equations to obtain 

B = i + f ^±|) 5 (5.25) 



a V« + L 
5 = A/2a (5.26) 

where the second equation is required since the right-hand-side of eqn (5.24) must 
be finite as V — > oo. Eliminating A gives 

B + f ^— -) B = (5.27) 



,1 + a 

and this leads to the general solution 

C 



A = A OQ + Cz-V- a M 1+a \ B = B^ + (j^j (5.28) 
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where A^, and C are integration constants. Eqns (2.17) and (2.18) give 
another relationship between these constants: 

-(i+aMc/a _ ( l + Ga + a 2 \ (1 + ag 6Boo . 2Aoo/a | 



C 



g2Boo-^Tu;jioo/« I ; I __|_ 



4a / 4a 



(5.29) 

where we have omitted a prefactor which depends on a. This shows that the 
asymptotically static solutions are described by two parameters for a given equa- 
tion of state. 

It should be stressed that the description "asymptotically static" in this con- 
text is rather misleading. This is because eqns (2.20) and (5.28) imply 

Vr « - 9 (1 ~ Q V ( 5 - 3 °) 
2a(l + a) 

at large z, so only the 1-parameter family of solutions with C = are asymp- 
totically static in the sense that the fluid is not moving with respect to the 
spheres of constant R. This agrees with the description of Foglizzo & Henrik- 
sen (1993), who term such solutions "symmetric". The exact static solution has 
= = C = 0, so the limit of the asymptotically static solutions need not 
be the exact static solution itself. We will describe the more general solutions 
with C 7^ as asymptotically "quasi-static" since they still have S and x tending 
to at infinity. 

At large values of z, eqn (2.20) implies that the energy function in these 
solutions is given by 



E = (1 + a) 2 e 6B 00 -2A 00 /a 

2(1 + 6a + a 2 ) 



l _ C(3 - a) ^_ (1 _ Q)/(1+Q) 
a(l + a) 



i. (5.31) 



where we have used eqns (5.28). The analogue of the parameter E in the dust 
solutions is therefore 

Eoo ~ e 6Bo °- 2Aoo/Q . (5.32) 

In the exact static case, the asymptotic value of E is just —M, where M is given 
by eqn (3.29). The analogue of the parameter D is more complicated but it is 
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related to the value of Vr at infinity. If one considers solutions in which E — > 
as z — > oo and then defines -D by analogy with eqn (4.19), eqns (5.30) and (5.31) 
imply 

D ~ e 3(l-a)Aoc/(2a)_ (533) 

It would be more natural to identify D with the value of z at the singularity (viz. 
Zoo = 1/D). However, this value can only be determined numerically and cannot 
be expressed in terms of the asymptotic parameters explicitly. 

To find asymptotically quasi-static solutions at small values of z, one seeks 
solutions with V = and finite values of A and B at z — 0. However, this 
requires A — B — 0, which from eqn (5.24) implies 

AB a = (^—^-) A . (5.34) 



a 

It is easy to see that this condition is incompatible with eqns (2.17) and (2.18), 
so there are no asymptotically quasi-static solutions at the origin (only the exact 
static solution itself). If instead we seek solutions in which A and B are finite 
and non-zero at z — 0, so that A and B diverge logarithmically, then eqn (5.24) 
requires 

A-fe). (,3 5 ) 

Substituting this into eqn (5.23) gives 

3(1 + a)B 2 + (5 + 3a)B + 2 = 0, (5.36) 
which has the two roots 

B = — ^ — - or B = -1. (5.37) 

3(1 + a) V ; 

However, these roots just correspond to solutions which are asymptotically Fried- 
mann or asymptotically KS at z — 0. We analysed such solutions in the previous 
sections, so the asymptotic quasi-static solutions do not yield any new behaviour 
at the origin. 
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The qualitative features of the asymptotically quasi-static solutions in the 
z > regime are indicated in Figures (6) for a = 1/3, Figures (7) for a < 1/3 
and Figures (8) for a > 1/3. As in the asymptotically Friedmann case, there 
are black hole solutions, a continuum of underdense solutions and discrete bands 
of overdense solutions (although the oscillations associated with the overdense 
bands are not exhibited in the figures). In each case, sample values of Eoo are 
shown and this illustrates the significance of the critical values E*, E crit and E stat . 

We have already emphasized that the z < solutions can be obtained from 
the z > ones by reflection about the origin. As in the dust case, the full family 
of solutions involve both regions. Indeed one has a figure precisely analagous to 
Figure (4), except that one now has a sonic point. Thus the collapsing solutions 
start in the z < region (following the reflection of the curves which go through 
the origin) and then jump to the z > region (following the curves which ap- 
proach a singularity). Likewise, the expanding solutions start off following the 
reflections of the curves which approach the singularity and then jump to the 
curves which pass through the origin. The fact that the solutions span both neg- 
ative and positive z has the important consequence that there may be more than 
one sonic point. The collapsing solutions have one sonic point in the z < regime 
but may have another two in the z > regime if the minimum of V is less than 
y/a. Likewise, the expanding solutions may have two sonic points in the z < 
regime if the maximum of V exceeds — ^fa and then another one in the z > 
regime. On the other hand, it is unlikely that solutions with more than one sonic 
point would be regular, so one would expect physical solutions to exist only for 
values of the parameters such that the maximum of \V\ exceeds \fa. 

A more extensive and quantitative analysis of the asymptotically quasi-static 
solutions, including a discussion of their regularity at the sonic point, is given 
elsewhere (Carr et al. 1998). In fact, some of the 2-parameter family of simi- 
larity solutions with pressure have already been studied numerically by Foglizzo 
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& Henriksen (1993), although they only focus on the collapsing solutions. [The 
relationship between their variables and ours is given in Appendix B.] They con- 
firm many of the qualitative features described above. In particular, they show 
that the solutions are described by two parameters at large \z\ and by one pa- 
rameter at small \z\ and they find the expected behaviour at the sonic point. In 
their phase space analysis, the orbits corresponding to the overdense solutions 
converge on and then spiral around the static solution for a while before head- 
ing to the origin. This corresponds to the oscillations found by Carr & Yahil 
(1990) and Ori & Piran (1990), with the number of oscillations identifying the 
overdensity band. Foglizzo & Henriksen (1993) confirm that the solutions with 
Vmin < 1 exhibit naked singularities. Indeed, the static attractor is closely related 
to the critical self-similar solutions found in the collapse calculations of Evans & 
Coleman (1995) and Maison (1996). This is discussed in more detail by Carr, 
Henriksen & Levy (1998). 

6 Discussion 

In this paper we have analysed the complete family of spherically symmetric 
similarity solutions for a perfect fluid with equation of state p = afi. The form 
of V(z) in these solutions is shown explicitly in Figures (6) to (8), although 
it should be emphasized that the figures are only qualitative and we have not 
included the (subsonic) solutions which are irregular at the sonic point. The key 
steps underlying our analysis are: (1) a delineation of the possible asymptotic 
forms at large and small distances from the origin; (2) an elucidation of the link 
between the z > and z < solutions; and (3) an explicit derivation of the dust 
solutions (which can be expressed analytically). 

In claiming that our classification is "complete", it should be emphasized 
that our considerations have been confined to similarity solutions of the first kind 
(i.e., homothetic solutions in which the similarity variable is z = r/t). However, it 
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may be possible to extend this work to the classification of similarity solutions of 
the "second" kind. For example, in spherically symmetric perfect fluid solutions 
which possess kinematic self-similarity, the similarity variable is of the form z = 
r/t a , where the a depends on some dimensional constant which contains a scale 
(Carter & Henriksen 1989). There is evidence that such solutions asymptote 
towards exact solutions that admit a homothetic vector (Benoit & Coley 1998), 
so the asymptotic analysis in this paper may be of rather more general application 
than is at first apparent; i.e. the asymptotic behaviour of all similarity solutions 
(not only those of the first kind) may be determined by the solutions described 
in this paper. Of course, the behaviour at finite values of the similarity variable, 
including for example the behaviour at sonic points and horizons, may be quite 
different. 

We have also confined attention to perfect fluids with a barotropic equation 
of state (necessarily of the form p = a/i) and so our analysis does not cover 
more general perfect fluids or anisotropic fluids, even though these may be of 
physical interest. In particular, a two-perfect-fluid model, in which each compo- 
nent is necessarily comoving and has an equation of state of the form p { = a^i 
(i = 1,2), is formally equivalent to a single perfect fluid that does not have an 
equation of state. It is therefore plausible that perfect fluid models for which pj ' \i 
is asymptotically constant may have the same asymptotic behaviour as the self- 
similar solutions studied in this paper. This is indeed the case for the two-fluid 
models in which each component separately satisfies the conservation equations 
(see Appendix C). 

We have not here considered the stiff case (a — 1), in which the speed of sound 
is equal to the speed of light. Since a = 1 is a bifurcation value, there can be 
significant changes in the qualitative behaviour from the a < 1 case. Therefore a 
discussion of stiff perfect fluid solutions may be important in understanding the 
dynamics of the complete class of similarity solutions. A partial analysis has been 
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made by Lin et al. (1976) and Bicknell & Henriksen (1978b). If a — 1, eqn (2.19) 
implies that V has no explicit dependence on z and, when V ^ 1, eqns (2.15) 
and (2.16) yield a single, second-order autonomous ODE for S. This equation 
can be better studied using different mathematical techniques to those employed 
in this paper (cf. Carr et al. 1998). In this context, it should be emphasized 
that our analysis does not cover the case in which the source is a massless scalar 
field since (if there is no scalar potential) this is formally equivalent to a stiff fluid 
whenever the gradient of the scalar field is timelike. The relevance of self-similar 
solutions to the occurrence of critical phenomena in scalar field collapse has been 
studied by Choptuik (1993), Abrahams & Evans 1993, Brady (1994), Gundlach 
(1995), Koike et al. (1995) and Frolov (1997). When a = -1, the perfect fluid 
source is equivalent to a cosmo logical constant. In this scale is introduced 

and so there are no self-similar solutions of the first kind. However, spherically 
symmetric self-similar solutions of the more general kind are still possible in this 
case (Henriksen, Emslie & Wesson 1983). 

Two special cases which have not been considered here are those in which 
the homothetic vector is either parallel or orthogonal to the fluid velocity. These 
cases are not covered by the analysis of Section 2. However, it can be shown 
that all perfect fluid spacetimes (not only spherically symmetric ones) admitting 
a homothetic vector parallel to the velocity vector are necessarily Friedmann 
(Coley 1991). In addition, Ponce de Leon (1993) has claimed that all spherically 
symmetric spacetimes which admit a homothetic vector orthogonal to the velocity 
vector have a singular metric (cf. Bogoyavlensky 1985). 

Acknowledgments 

We thank Martin Goliath, Dick Henriksen, Ulf Nilsson and Claes Uggla for 
useful discussions and Andrew Whinnett for help with the numerical work. An 
earlier version of this paper omitted the asymptotically Minkowski solutions and 



51 



we extended our analyis to cover these only after their existence was ascertained 
by Goliath, Nilsson and Uggla numerically. BJC is grateful to the Department 
of Mathematics and Statistics at Dalhousie University for hospitality received 
during this work. 



52 



References 

A. M. Abrahams and C. R. Evans, 1993, Phys. Rev. Lett. 70, 2980. 
G. I. Barenblatt and Ya B. Zeldovich, 1972, Ann. Rev. Fluid Mech. 4, 285. 
P.M. Benoit and A. A. Coley, 1998, Class. Quantum Grav. 15, 2397. 
E. Bertschinger, 1985, Ap. J. 268, 17. 

G. V. Bicknell and R. N. Henriksen, 1978a, Ap. J. 219, 1043. 

G. V. Bicknell and R. N. Henriksen, 1978b, Ap. J. 225, 237. 

G. V. Bicknell and R. N. Henriksen, 1979, Ap. J. 232, 670. 

O. I. Bogoyavlenski, 1977, Sov. Phys. JETP 46, 634. 

O. I. Bogoyavlenski, 1985, Methods in the Qualitative Theory of Dynamical 
Systems in Astrophysics and Gas Dynamics ( Springer- Verlag). 

P. R. Brady, 1994, Class. Quant. Grav. 11, 1255. 

A. H. Cahill and M. E. Taub, 1971, Comm. Math. Phys. 21, 1. 

B. J.Carr, 1995, preprint prepared for but omitted from The Origin of Struc- 
ture in the Universe, ed. E.Gunzig and P.Nardone (Kluwer). 

B. J. Carr, 1998, in Proceedings of the Seventh Canadian Conference on 
General Relativity and Relativistic Astrophysics, ed. D. Hobill. 

B. J. Carr and A. A. Coley, 1998a, Class. Quant. Grav., in press. 

B. J. Carr and A. A. Coley, 1998b, Class. Quant. Grav., in press. 

B. J. Carr and S. W. Hawking, 1974, MNRAS 168, 399. 



53 



B. J. Carr and A. Koutras, 1993, Ap. J. 405, 34 (CK). 

B. J. Carr and A. Yahil, 1990, Ap. J. 360, 330 (CY). 

B. J. Carr and A. Whinnett, 1998, preprint. 

B. J. Carr, R. N. Henriksen and M. Levy, 1998, preprint. 

B. J. Carr, A. A. Coley, M. Goliath, U.S. Nilsson and C. Uggla, 1998, 
preprint. 

B. Carter and R. N. Henriksen, 1989, Ann. Phys. Supp. 14, 47. 
M. W. Choptuik, 1993, Phys. Rev. Lett. 70, 9. 

M. W. Choptuik, 1994, in Deterministic Chaos in General Relativity, ed. 
D. Hobill et al. (Plenum, New York). 

A. A. Coley, 1991, Class. Quant. Grav. 8, 955. 

A. A. Coley, 1997. in Proceedings of the Sixth Canadian Conference on 
General Relativity and Relativistic Astrophysics, The Fields Institute Com- 
munications Series (AMS), Volume 15, page 19, eds. S. P. Braham, J. D. 
Gegenberg and R. J. McKellar (Providence, R.I.). 

A. A. Coley and B. O. J. Tupper, 1989, J. Math. Phys. 30, 2616. 

C. B. Collins, 1977, J. Math. Phys. 18, 2116. 

D. M. Eardley, 1974, Comm. Math. Phys. 37, 287. 

C. R. Evans and J. S. Coleman, 1994, Phys. Rev. Lett. 72, 1782. 
T. Foglizzo and R. N. Henriksen, 1992, Phys. Rev. D. 48, 4645 (FH). 
A. V. Frolov, 1997, preprint. 

54 



M. Goliath, U. S. Nilsson and C. Uggla, 1998a, Class. Quant. Grav. 15, 
167. 

M. Goliath, U. S. Nilsson and C. Uggla, 1998b, Class. Quant. Grav. 15, 
2841. 

C. Gundlach, 1995, Phys. Rev. Lett. 75, 3214. 

R. S. Hamade, J. H. Home and J. M. Stewart, 1996, Class. Quant. Grav. 
13, 2241. 

R. N. Henriksen and K. Patel, 1991, Phys. Red. D. 42, 1068. 

R. N. Henriksen and P. S. Wesson, 1978, Astrophys. Space Sci. 53, 429. 

S. Ikeuchi, K. Tomisaka and J. Ostriker, 1983, Ap. J. 265, 583. 

P. S. Joshi and I. H. Dwivedi, 1993, Phys. Rev. D. 47, 5357. 

R. Kantowski and R. Sachs, 1966, J. Math. Phys. 7, 443. 

T. Koike, T. Hara and S. Adachi, 1995, Phys. Rev. Lett. 74, 5170. 

K. Lake, 1992, Phys. Rev. Lett. 68, 3129. 

R. B. Larson, 1969, MNRAS 145, 271. 

D. N. C. Lin, B. J. Carr and S. D. M. Fall, 1978, MNRAS 177, 151. 
D. Maison, 1996, Phys. Lett. B. 366, 82. 

B. D. Miller, 1976, Ap. J. 208, 275. 

C. W. Misner and H. S. Zapolsky, 1964, Phys. Rev. Lett. 12, 635. 
U. S. Nilsson and C. Uggla, 1997, Class. Quantum Grav. 14, 1965. 



55 



A. Ori and T. Piran, 1987, Phys. Rev. Lett. 59, 2137. 

A. Ori and T. Piran, 1988, Mon. Not. R. Astron. Soc. 234, 821. 

A. Ori and T. Piran, 1990, Phys. Rev. D. 42, 1068. 

M. V. Penston, 1969, MNRAS 144, 449. 

J. Ponce de Leon, 1988, J. Math. Phys. 29, 2479. 

J. Ponce de Leon, 1990, J. Math. Phys. 31, 371. 

J. Schwartz, J. P. Ostriker and A. Yahil, 1975, Ap. J. 202, 1. 

L. I. Sedov, 1967, Similarity and Dimensional Methods in Mechanics (New 
York, Academic). 

J. Wainwright and G. F. R. Ellis, 1997, Dynamical Systems in Cosmology 
(Cambridge University Press, Cambridge). 

P. Wesson, 1986, Phys. Rev. D 34, 3925. 

P. S. Wesson, 1989, Ap. J. 336, 58. 

A. Whitworth and D. Summers, 1985, MNRAS, 214, 1. 

Appendix A. 

In this paper we have used "comoving" coordinates, since this approach is best 
suited to studying the solutions explicitly. However, it should be stressed that our 
work is complemented by the analysis of Bogoyavlensky (1985) and Goliath et 
al. (1997, 1998) using "homothetic" coordinates and that of Ori & Piran (1990) 
and Maison (1995) using Schwarzschild coordinates. In this Appendix, we discuss 
these other approaches in more detail. 



56 



In the nomothetic approach, the coordinates are adapted to the homothetic 
vector and this yields results which complement and, in some cases, provide more 
rigorous demonstrations of the conclusions reached in this paper. However, in the 
homothetic approach, spacetime must be covered by several coordinate patches, 
one in which the homothetic vector is spacelike and one in which it is timelike. 
These regions must then be joined by a surface in which the homothetic vector is 
null and this surface is associated with important physics. Bogoyavlensky (1985) 
studied the spacelike and timelike cases simultaneously (with the metric being 
written in "conformally static" form) and continuously matched the two regions 
to obtain the behaviour of solutions crossing the null surface. However, it should 
be noted that Bogoyavlensky changed comoving coordinates explicitly to describe 
the physics of the associated solutions. 

Recently Goliath et al. (1998a, 1998b) have reinvestigated both the spa- 
tially and temporally self-similar cases. The timelike region contains the more 
interesting physics (eg. shocks and sound-waves). They introduce dimensionless 
variables, so that the number of equations in the coupled system of autonomous 
differential equations is reduced, with the resulting reduced phase-space being 
compact and regular. In this way the similarities with the equations governing 
hypersurface orthogonal models, and in particular spatially homogeneous models 
(Wainwright & Ellis 1997, Nilsson & Uggla 1997), can be exploited. In their 
approach, all equilibrium points are hyperbolic, in contrast to the earlier work 
in which Bogoyavlensky (1985) used non-compact variables (which resulted in 
parts of phase-space being "crushed"). One drawback with this approach is that 
spacetime must be covered by more than one patch and the solutions must be 
matched between the spatially and temporally self-similar regimes in order to 
determine their physical properties. 

The Schwarzschild approach is better suited to studying the causal structure 
of the self-similar solutions. This is because, in order to obtain physically rea- 
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sonable models, spacetimes are often required to be asymptotically flat. Since 
asymptotically flat spacetimes are not self-similar, one therefore needs to match 
a self-similar interior region to an non-self-similar exterior region and this is usu- 
ally taken to be Schwarzschild. In particular, Schwarzschild coordinates are most 
suitable for solving the equations of motion for (radial) null geodesies, as required 
in studying the global structure of the solution. Consequently it was used by Ori 
& Piran (1990) since one of their primary goals was to study naked singularities 
and test the cosmic censorship hypothesis. However, the Schwarzschild coordi- 
nates break down at t — 0. 

Appendix B 

The precise transformations between the various coordinate systems used to 
study self-similar spherically symmetric perfect fluid models are given explictly 
in Bogoyavlenski (1985; see Section 3 of Chapter IV). The coordinate transfor- 
mations between the comoving and Schwarzschild systems, both of which are 
employed by Ori & Piran (1990), are given explicitly in their paper. The trans- 
formations between the homothetic and Schwarzschild coordinates and between 
the homothetic and comoving coordinates are given explicitly in Goliath et al. 
(1998a; Appendix B), where the relationship between their variables and those 
of Ori & Piran (1990), Maison (1996) and Foglizzo & Henriksen (1992) are also 
given. The relationsip between their variables and those used by Bogoyavensky 
(1985) are given by Goliath et al. (1998). 

Here we explicitly demonstrate the relationship between the variables used 
in this paper and those used in Foglizzo & Henriksen (1992; FH). The main 
functions used in FH are the three functions (N, ~p, V 2 ), defined by eqns (FH3) 
- (FH5), which depend on the similarity variable £ = z~ x . The remaining self- 
similar functions can then be written in terms of these [see eqns (FH6) - (FH8)]. 
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Their function V is identical to ours. Using eqns (2.12) and (FH3) we find that 

N(z) = ai :r (a - 1)A V {Q ~ 1)/{1+Q) . (Bl) 



Using eqns (2.8), (2.12), (2.17) and (FH4), we obtain 

S 



fi(z) = 3 



l + (l + a) s 



(B2) 



Conversely, x and S can be defined explicitly in terms of (JV, /x, V 2 ) through eqn 
(FH8), 

S 2 = a 2 z^/^[N\V\}-\ (£3) 

and eqn (Bl). 

The differential equations governing the evolution of (JV, Jl, V 2 ) are given by 
eqns (FH12) - (FH14); these equations constitute an autonomous system of ODEs 
in terms of the variable In^ = —Inz. Eqns (FH12) and (FH13) are equivalent 
to our eqns (2.15) and (2.16). The first integral of the governing ODEs is given 
by eqn (FH10) and is equivalent to our eqns (2.17) and (2.18). Eqn (FH14), 
which governs the evolution of V 2 , is obtained by differentiating V 2 , defined by 
eqn (2.19), and using the first integral. Consequently, the evolution eqn (FH14) 
replaces eqns (2.17) and (2.18). 

FH then regularize their system of eqns by introducing a new independent 
variable, r, defined by eqn (FH15), which is equivalent to 

dlnz o . „ JN 

—— = -1 + ay- 2 . (754) 
dr 

This divides phase-space into two disconnected components. Although the re- 
sulting system of ODEs is autonomous, the system is not regular at £ = despite 
the fact that this point does not correspond to a physical singularity (i.e., it arises 
due to a coordinate problem at £ = 0). FH then introduce new functions and 
coordinates so that solutions are completely regular at (t — 0, r > 0). However, 
the resulting ODEs are no longer autonomous after this transformation. 
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Appendix C. 

The expansion of the comoving fluid velocity congruence, 6 = u a - a , is given by 

rO = ze- u Q (CI) 

where 

e(z) = -^-(\ + 2S). (C2) 

dz 

When p = afj,, the conservation equations then yield 

dW 



dz 



= -(l + a)WQ. (C3) 



If we consider two comoving perfect fluids as the source of the gravitational field, 
each of which satisfies 

Pi = Oifii, Wi = fiiR 2 {a = 1, 2) (C4) 

(e.g., a mixture of dust and radiation with ct\ = 1/3 and a-i — 0), then the source 
is equivalent to a single perfect fluid with 

A* = A*i + A*2, P = Pi + P2 = oci fix + a 2 ^ 2 (C5) 

although this does not admit an equation of state. 

Suppose the two perfect fluids are non-interacting, with each separately sat- 
isfying the conservation equation (C3). Then 

dW, 
dz 

We define a new variable, x> by 



-(l+a t )W t Q- (C6) 



/il - /i 2 Wi - 

* = xW = = mTrn- {C1) 

where — 1 < x < 1. From eqn (C6) we derive the evolution equation for x' 

^ = ^(«i-« 2 ){l-X 2 } (C8) 
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where r is defined by 

for regions in which 9 (and hence the expansion 6) is non-zero. 

Eqn (C8) is a decoupled autonomous equation for x- It has equilibrium points 
X — il, and hence all solutions asymptote to % = ±1 in regions for which the 
expansion does not become zero, x = +1 corresponds to fi2 = and x = ~ 1 
corresponds to fii = 0; i.e., the solutions of these self-similar two-fluid models 
asymptote towards the exact asymptotes of the single perfect fluid solutions. 
Asymptotically 

p/fi-^ai, (CIO) 

and which value of ctj is picked out (i.e., which of the two single fluids govern 
the dynamics asymptotically) depends on the signs of (cti — a 2 ) and 6 and on 
whether \z\ — > or \z\ — > 00. 

Figures 

FIGURE (1). This shows the form of V(z) for the exact Friedmann (F), static 
(S) and Kantowski-Sachs (KS) solutions for (a) the general a case and (b) the 
a = 1/3 case. Also shown are the sonic lines |V| = \ j\foL (broken), the range of 
values of z in which curves can be regular at the sonic point (z\ < z < z 2 and 
z > z 3 ) and (for a — 1/3) the curves corresponding to M = (solid). 

FIGURE (2). This shows the form of the scale factor S(z), the velocity 
function V(z) and the density function fit 2 (z) for the asymptotically Friedmann 
dust models. These are described by a single parameter E where E = in the 
exact Friedmann case: the z > solutions are overdense and collapse to black 
holes for E < (with an event horizon for E > E* since V min < 1) but they are 
underdense and expand forever for E > 0. The z < solutions are just the time 
reverse of these. The broken curve corresponds to a solution with negative mass 
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and is probably unphysical. 

FIGURE (3). This shows the form of the scale factor S(z), the velocity 
function V(z) and the density function nt 2 (z) for dust models with E — and 
different values of D. These solutions necessarily span both positive and negative 
values of z. For D > they represent expanding (solid) or collapsing (broken) 
solutions and the latter contain a naked singularity (V m i n < 1) if D exceeds 
D*. The D < models (dotted) undergo shell- crossing before encountering the 
singularity and are probably unphysical. 

FIGURE (4). This shows the form of the scale factor S(z), the velocity 
function V(z) and the density function fit 2 (z) for the asymptotically quasi-static 
dust solutions. These are described by two parameters (D and E). For E > 
they resemble the curves in Figure (3), with both expanding (solid) and collapsing 
(broken) solutions. The collapse singularity is naked (V min < 1) if E exceeds a 
value E*(D). For E < there are also solutions which recollapse to black holes 
or emerge from white holes (broken), as in the asymptotically Friedmann case. 
The dotted curve corresponds to an (unphysical) negative mass solution. 

FIGURE (5). This shows the form of the asymptotically Friedmann solutions 
for a radiation equation of state (a = 1/3), with particular emphasis on the 
behaviour at the sonic point. Solutions which are regular (irregular) at the sonic 
point are shown by solid (dotted) lines, while black hole solutions (with no sonic 
point) are shown by broken lines. There is a 1-parameter continuum of underdense 
solutions which are regular at the sonic point but the overdense solutions lie in 
discrete bands (just the first of which is shown) and are characterized by the 
number of oscillations they exhibit. 

FIGURE (6). This shows the form of the velocity function V(z) for the full 
family of spherically symmetric similarity solutions with a = 1/3. The exact 
Friedmann, Kantowski-Sachs and static solutions are indicated by the bold lines. 
Also shown are the asymptotically Friedmann solutions (for different values of E) 
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and the asymptotically quasi-static solutions (for different values of E and D). 
The asymptotically quasi-static solutions contain a naked singularity when the 
minimum of V is below 1. The negative V region is occupied by the asymptoti- 
cally Kankowski-Sachs solutions, though these may not be physical since the mass 
is negative; solutions which are irregular at the sonic point are shown by broken 
lines. Also shown are some of the solutions which are asymptotically Minkowski 
at finite and infinite values of z. 

FIGURE (7). This shows the form of the velocity function V{z) for the full 
family of spherically symmetric similarity solutions with a < 1/3. 

FIGURE (8). This shows the form of the velocity function V(z) for the full 
family of spherically symmetric similarity solutions with a > 1/3. 
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